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AN INTRODUCTION TO THE LAWS OF AIR RESISTANCE OF AEROFOILS. 



By George db Bothezat. 



PREAMBLE. 

It is only very slowly, tlu-ough the centuries, that the notion of the resistance of a fluid 
to the motion of a solid body has been developed. This notion is intimately associated with 
tlie concepts which we gain from mechanical ])henomena. In the aurora of the first gleams 
which pierced the darkness of the human mind in the domain of the concepts of motion, fluid 
resistance was not differentiated from motion. Thus Aristotle ' considered in principle — not 
willing to admit the possil)ility of a vacuum — that the resistance of a fluid was inseparable 
from the phenomena of motion. It is this point of view which paralyzes, so to say, completely 
his attempts to form a conception of the phenomena of motion, the exposition of which by 
him was, it must be adfled, very hazy. Through antiquity to the Middle Ages, dynamical 
phenomena were dawning, but with a very confused misundci-standing. Leonardo de Vinci 
seems to have thought much about the motion of bodies under terrestial conditions. It is 
without any doubt that he made numerous and remarkable attempts at mechanical flight. 
But, in his dynamical concepts, he does not seem to have clearly separated the phenomena of 
motion from the phenomena of the resistance of fluids. Thus he used the confused conception 
of the impetus which ought to be communicated to a body when the same is set in motion, 
and which ought to dissipate itself progressively to cause the body to stoj). But by the use of 
the conception of dissipation of the impetus, he even arrived at the happy conclusion of the 
impossibihty of perpetuum mobile. It is Gahleo^ who finally has a full conception of the material 
nature of the gas&s and of the influen(^e the same have on the motion of bodies — an influence 
which he knew to decrease with the velocity. This is why Galileo, in his celebrated experi- 
ments on falling bodies, recognized tlie necessity of making them at low velocities. Low veloci- 
ties first made possible the quantitative observations, and secondly diminished all the resist- 
ance, for the decreasing of which all possible measures were taken. So it is that Galileo first 
came to the modern <'onception of dynamical phenomena. 

To disengage the law of the motion of bodies, considering the latter as moving in vacuo 
and without any kind of resistance, and to look on all other effects, such as friction or medium- 
resistance, as additional effects, this was the conception which allowed the establishTnent of 
dynamics. ITiis conceptional sorting of questions in the complex problem of motion must be 
considered as one of the greatest scientific conquests. On our planet the motion of bodies 
always takes place in a fluid. The phenomena of motion, taken as a whole, is so complex that 
it is inextricable for the human mind. A very large conceptional effort had to be developed to 
rise to the abstraction of the phenomena of motion, cleared fi'om the influence of the immediate 
medium. But once this big step made, we have the magnificent picture of the powerful dynami- 
cal laws, of which we have seen the development; the questions of friction and fluid resistance 
being considered as special separate questions, whose complexity is enough to make them subjects 

■ Aristotle. ■' Physique," Livro VH, Chap. VI. French translation Ijy ISarthelniy Saint-Hilaire. 

< Galileo. "Dlscorsi et dimostrazioni matcmaticbc intomo a due nuo\c scienie attenanti alia mccanica ct 1 movimciiti looaU." Bologna 
edition 1655. 
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of special branches of science. We understand now why all the first dynairdcal investigations, 
until recent times, were almost exclusively made on problems wherein the influence of the 
medium could be neglected. It is for this last reason that the advance of celestial mechanics 
was developed long before we began to understand the laws of the motion of bodies under terrestial 
conditions. It is undoubtedly true that it is only with the birth of a\aation that j^rofound 
studies of the motion of bodies in real fluids were started and that the light began to penetrate 
through the complex and delicate phenomena of fluid resistance, which phenomena have, for 
a long time, veiled from our eyes the laws of dynamics but which have now given us the con- 
quest of the aerial ocean. 

Actually we are only taking the first steps in the conception of the problems of fluid resistance. 
The fornier status of these questions consisted more in the comparison of fluids to some mechani- 
cal system, more or less similar to fluids, than in the studj" of the real fluids with their real prop- 
erties. Thus, Newton likened fluids to a system of elastic particles whose impact on the soUd 
body produced the fluid resistance. Euler ' in his research on fluid resistance, likened fluids 
to a continuous homogeneous frictionless medium and calculated the fluid resistance ])y aid of 
the general equations which he built up for that kind of medium. He was brought to the con- 
clusions, very far from reality, that a body moving in a fluid meets no resistance to its motion. 
This conclusion is a consequence of the assumption of a continuous and noncyclic flow around 
the body. Recently Kutta has shown that, in the general case of the continuous flow of a 
perfect fluid around the body, the circulation around tlic contour embracing tlie body can have 
a finite value, and in such a case the fluid resistance has a finite value but is perpendicular to 
the general stream velocity. Thus, in a perfect fluid only the power corresponding to the 
resultant pressure on the surface of the body is necessarily equal to zero, but the resultant 
pressure can have a finite value. We will later consider Kutta's conceptions. Helniholtz is 
the first to have made a serious attempt to bring the foundation of hydrodynamics into more 
close agreement with reality; and his work in that sense is of gi-eat iniportance. He showed 
the necessity for the consideration of vortex nrotion and indicated the possibihty of the forma- 
tion of surfaces of discontinuity in fluid motion. 

Tliis last idea of surfaces of discontinuity was used by KirchhofF ' and by Lord Rayleigh' 
for the calculation of fluid resistance in some simple cases, which method was recently 
largely developed by G. Greenhill * H. L«vy ° and others. The flow which in reality is estab- 
lished seems only rarely to be of the kind assumed by Kirchhoff and Lord Rayleigh, so that in 
general the experimentallj' measured fluid resistance does not corresjiond to that calculated by 
the Kirchhoff and Lord Rayleigh method, a fact to which already WiUiam Thomson ' (Lord 
Kelvin) has drawn attention. The way in which viscosity has been considered until now 
does not give a satisfactory solution of the problem of fluid resistance, either. The calculation 
of fluid resistance by the equations of motion of a viscous fluid in the final form given them by 
Stokes ' seems to agree with e.xperiment only for very small velocities. It is only the develop- 
ment of aviation that has given a new powerful impulse to aerodynamics, and has brought with 
it the necessity of a conception of fluid resistance closer to reality. Many c^uite new ideas and 
concepts have thus been progressively developed. 

In 1902 W. M. Kutta ' formulated, first for a particular case and soon after generalized for 
the general case, an important theorem which gives the relation between the fluid resistance 
and the flow around a l)ody which encounters that resistance. This theorem was established 
by its author for the case of perfect fluids. In that case, this theorem tells us that the lift of 

> Euler. "Prlncipes generaux du mouvemcnt dcs nuidcs." Eistolre dc I'Acad. de Berlin 1755. "Cc principiis motus Fluldorum. " Novl 
Comm. jicad. rctrograd 1759. 

J Kirehhotl. " Vorlcsimgen abcr Mcchanik," 1&97. 22st Vorlcsung. 
» Lord Rayleigh. "Scientific papers," I. 

• G. nrccnhill. "Stream lines piist a plane barrier, and oftbe discontinuity arising at an edge." Beport 19, Advisory Committee lor Aeronau- 
tics, 1912. 

• H. Levy. "Discontinuous motion past a Curved Boundary," Proc. Koy. Soc. 1916. 

• Lord Kelvin. "Math, and I'hys. papers." IV, p. 215. 
' Stokes, G. O. " Math, and Phys. papers. " I. p. "8. 

' W. M. Kutta. "lUustrirte Aeronautische Mittailungen," 1902. "Sitzungsberichte der KonigUchen Bayeriscben Akademie der Wissen- 
schaltan." Milnchen, 1910 and 1911. 
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the body is equal to the product of the density, velocity of the stream at infinity, and circu- 
lation around a contour surrounding the body, but for the drag it gives the value zero. To cal- 
culate the lift by this theorem, it is necessary to know the flow around the body. Kutta 
assumes that this flow is a compound forward and cyclic motion of the fluid; but this is only 
an assumption, without sufficient foundation. The Kutta theorem, understood as just stated, 
was applied with many developments to numerous cases by Joukowski ' and Tchapliguine. 
The results of all these calculations do not fully agree with experunent. I have siibmitted the 
Kutta theorem to critical examination ^ and have showed that this theorem must not be under- 
stood as giving the solution of the problem of fluid resistance, because it leaves open the ques- 
tion of the flow around the body and only gives the relation between flow and fluid resistance. 

In recent years, Karman ' has called attention to the fact that the flow around a body 
having a rectilinear and uniform motion of translation in a fluid very often consists of a system 
of vortices whicli are formed behind the bod3% and has shown the relation which must exist 
between the momentum of these vortices and the fluid resistance.* 

The works of Karman are probably the first to indicate the necessity of the determination 
of the type of flow wliich in reality takes place around a body in order to bo able to calculate 
ts fluid resistance, a question to which not enough attention was paid before. And what is 
particularly important, the type of flow which most generally establishes itself is not necessarily 
one of the types which were presupposed by all the foregoing theories. The flow around a 
body immersed in a fluid is not necessarily continuous as it was supposed by Euler; it is not 
generally characterized by a system of surfaces of discontinuity either, as was assumed by 
Kirchhofl' and Lord Rayloigh, wliich surfaces of discontinuity must be considered as almost 
unstable, the viscosity disturbing them; but more often the flow is characterized by a system of 
vortices as shown by Karman. 

Nevertheless, the systematical study of the different kinds of flow around solid bodies 
which are compatible with the general equation of hydrodynamics is of the highest value. It is 
of the geratest importance to disengage all the types of flow which are possible for fluids because 
only under such conditions can we reach the complete solution of the great fluid resistance 
problem. Generally speaking, all kinds of flow satisfying the equations of hydrodynamics are 
virtually possible under special conditions. Particular attention must, however, be paid to 
the question of finding out the exact conditions under which each kind of flow can take place. 
In many cases of flow of air or w"ater the types of flow characterized by vortices in quincunx 
seem to be most usually obtained. This is on account of the need of stability, and the con- 
ditions of energy dissipation inside those fluids. If we look over the historical development 
of hydrodynamics it is the progressive discovery of the properties of the different types of flow 
that we see before us. 

In the development of modern hydrodynamics the question of the conditions which fix 
the type of flow established under given conditions was left nearly without any examination. 
Exactly spoakuig, what did the classical hydrodynamics give us m order to detennine the flow 
in the case of steady motion ? Of the four equations of the motion of an incompressible fluid 
which forms the foundation of classical h3'drodynamics, three give the relation between the 
distribution of the velocities and the pressures ' — it is these which express the theorem of 
momentiun in its application to a fluid particle — and only one, the equation of continuity 
determines the flow. The question of finding the flow around a body as defined by the equa- 
tion of continuity is a problem of finding a function which verifies the Laplace equation and 
satisfies the boundary conditions. It must be remembered, however, that the equation of 
continuity is oidy a necessary condition for continuity and is not at all sufficient. As Helm- 
holtz has first remarked, the discontinuity of the tangential components of the velocity in 

' See "AerodjTiamique" by N. Joukowski. Paris, 1916. 
2 See Note I at the end of this pamphlet. 

'V. Karman. " N'aohrichtcn von der KOnigllchen Gesellschalt der Wissenschaften lu Gdttlngen." 1911. "Pbysikallsche ZeitscbiUt." 1912. 
See also the above-mentioned Aerodynamique, by N. Joukovrski and Note IV' at the end of this pamphlet. 
* In the following the conceptions of Karman will be extended to the aerofoil. 

' This Is particularly well seen when we use the equations of Uuid motioD in natural curvUineor coordinates. See Note II. 



10 



ANNUAL REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 



regard to some surfaces is compatible with the equation of hydrodynamics, so that when a flow, 
satisfying the equation of continuity is found, it must still be verified that such a flow is vir- 
tually possible. This is probably one of the most important questions in the problem of fluid 
resistance. I must finally add that in some cases the continuous flow of a fluid seems to be 
practically impossible. 

I will give an example. To follow more easily the motion of a fluid, let us divide its con- 
tinuous volume by a system of triorthogonal surfaces which accompany the fluid in its motion, 
so that we have no flow through these surfaces. Continuous motion wiU mean that each 
fluid element will always remain in contact with the 14 elements which are in touch with it 
at any moment; that all the elements coutamed in any closed surface moving with the fluid 
will always remain in it; that all the elements which are inside the fluid will never come on its 
surface; that all the elements which are on the boundary surface of the fluid will never come 
inside the fluid, etc.; so that the whole motion is consideretl only as a continuous deformation 
of the fluid medium without any alteration of the mutual grouping of the elements. If we now 
consider for example the flow of a viscous fluid running out of a pipe into a reservoir, considering, 
as generally admitted, the velocity of the fluid on the pipe walls equal to zero, and if we attempt 
to follow the deformation of a fluid element, we very easily see the impossibility of such a concep- 
tion. It is enough to remember that the elements all keeping close togetlier will be found in some 
cases making some hundreds of thousands of revolutions per second.' The admittance of conti- 
nuity in such conditions seems to be very difficult. In aU probability, the real motion must 
consist of a succession of continuous states of motion interrupted by discontinuous intervals. 

The following question can very natiu-ally arise: How did it happen that in the domain 
of rigid dynamics we at once reached so many results whicli stay in close agreement with the 
motion of real solid bodies, and tliat in many hydrodynamical problems we have not been till 
now able to secure satisfactory solutions. The fact lies in the nature of the question. In 
the historical evolution of mechanics the concept of a rigid body was first fully reached. The 
formation of this concept did not present an}' special difl^culties and its application to the 
analysis of an enormous number of problems of practical mechanics has shown at once all its 
power. The scientific world was already in the possession of a fully developed rigid dynamics, 
experimentally verified, when, in Euler's times, attention was brought to the general problem 
of fluid motion. When the concept of a perfect fluid was readied it was instinctively assimaed 
that this conception bore a relation to the real fluid quite as close as the conception of rigid 
bodies to a real solid body. It was with great astonishment that men recognized the disagree- 
ment which began to appear between tlie consequences of the hydrodynamical equations and 
the hydraulic experiments. For a long time the investigators in hydrodj^namics somewhat 
skeptically considered the disagreement between theory and practice, and did not pay much 

• In the case of a continuous motion of a fluid in a horizontal pipe, the axis of the pipe is an axis of summotry for the whole phenomenon. Using 
this axis, as the Z axis of a system of cylindrical coordinates, the equations of the Quid motion in these coordinates for our case wiU be 

bv, a«t), dp 

Where s is the density of the fluid and m is the viscosity constant. But in consequence of continuity and incompressibility 

5j!!_0and^?-.4-Const. 



where A is the pressure gradient along the Z axis, so tliat 
and 

For the vortex components wo have 
and tram the foregoing 
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In some of Poiseuille's experiments with water, A was of the order of one atmosphere s IC ; csO.Ol ; r~ 0. 01 
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attention to it; and I will allow myself to say that probably some of them believed more the 
conclusions of their equations rather than the experimental results. The great success of rigid 
dynamics in its origin is without any doul)t one of the principal reasons why there has been 
so much confidence in the concept of a perfect fluid and wliy in a certain period of the devel- 
opment of hydrodynamics this science has been brought to a very abstract development, 
more as a mathematical discipline than as a science of nature. But tlie demands of the mag- 
nificent conquest of the aerial ocean by the airplane has, I think, definitely brought the 
hydrodynamical science on the right way of one of the most important natural sciences. 

If we review the foregoing, we can now give the foUomng statement of the question of 
fluid resistance: 

To be able to calculate the fluid resistance of a body, we must first determine the type 
of flow which takes place around tlie body in the case considered. 

It appears that the conditions which hold at the surface of contact of fluid and solid con- 
stitute a special difliculty and hence that special conditions exist there. It may therefore 
seem that it is necessary first to make a special studj- of the problem of the floAv of the fluid 
in the immediate neighborliood of tlie body. As conditions of flow de])end upon tlie sliape 
of the surface of the body and the physical properties of the fluid, I think that the solution 
of this problem could be obtained only in an empirical-theoretical way; that is to say, to find 
out by what quantities, experimentally measured, we can fix tlie mutual relation between the 
surfaces of contact and the fluid flowing along them, so tliat these quantities once kno'WTi, 
the flow in the neighborhood of the body could be determined. It seems that only a thin 
layer of fluid is disturbed by the immediate influence of the surface of a body and that at a 
moderate distance from the body the influence of the body surface practically disajjpears. 

The conditions of flow in the portion of a fluid remote from any rigid body seem to be 
easier to imderstand than the conditions in the immediate neighborhood of a body. In the 
remote fluid portions we can have continuous motion, and so long as continuous motion takes 
place no vortices can appear within the fluid, and this independently of any assumption as to 
viscosity.^ The appearances of vortices can only come from the formation of surfaces of dis- 
continuity in the fluid. The mechanics of formation of the latter surfaces is very probably 
the following: 

A real fluid has to be considered as a fluid-elastic body (in opposition to the solid-elastic 
body), the stresses in which are fixed by the distribution of the velocity gradient. The fluid- 
elastic body can, without any doubt, move as a continuous whole only provided the stresses 
at all the points of the fluid have not reached a certain value. If some of these stresses exceed 
a certain magnitude, which must depend upon the properties of the fluid, the fluid may break 
at that point if tension stresses appear, or slip, if the stresses are shears. It is in this way 
that surfaces of discontinuity arise in a fluid. But the existence of them can be only a momen- 
tary phenomenon which is replaced by vortices, the surfaces of discontinuity being unstable 
in regard to viscosity. We thus see that the study of the problem of fluid resistance must 
consist first, of finding out the conditions imder which continuous motion of a fluid can take 
place around a body. The sj'stem of stresses in the fluid aroimd the body seems to be the 
criterion for that continuity. When the latter conditions are not satisfied, then we shall have 
to find out what systems of vortices can be compatible with the problem; then, afterwards, 
when the type of flow is exactly fixed, the fluid resistance can be calculated by the theorem 
of momentum. 

I Let lis cousider. for simplicity, a viscous fluid moving parallel to a plane. (For the notations, SCO Note II.) The moment olthe forces acting 
on two opposite sides of a fluid element will be 

^^^^ , AO^' 

that is, an uifinitcly small quantity of the fourth order. 

The moments of momentum of the element will he proportional to 

that is, an inlluitoly small quantity of the fifth ordi>r. Tlicrofore, if the forces acting on the surface of the element give rise to a moment, the 
element will necessarily lake an inDmte rotation, wliich would mean discontinuity. 
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I will here ask the reader to pay special attention to the following fact. The formation 
of special types of vortex systems behind a body moving in a real fluid is a direct consequence 
of the energy dissipation inside the fluid owing to viscosity. As will be shown in this pam- 
phlet, tlic work of the fluid resistance forces brings with it the necessity of the formation behind 
the body, in the limiting case, not only of the interior surfaces of discontinuity of the Kirchhofl- 
Lord Rayleigh theory, but also of a system of exterior surfaces of discontinuity remote from 
the body. All these surfaces of discontinuity are constituted by vortex sheets. But such 
vortex surfaces of discontinuity being imstable, thej- go over into stable vortex systems, the 
quincunx vortex sj'stem being the one most generally obtained, for tlie case of large aspect 
ratio. I call fundamental wave the vortex motion generated by the exterior vortex surface of 
discontinuity, and secondary wave the vortex motion generated by the interior vortex surface 
of discontinuity. We can now understand why at small velocities the flow around a body 
approaches more a continuous flow. At small velocities tlie work of the fluid resistance forces 
is small and is quickly dissipated inside the fluid. But at greater flow velocities the work of 
the fluid resistance forces can not be at once dissipated in the fluid, and a decrease of the kinetical 
and potential energy of the fluid is produced, which gives rise to an oscillator}' motion of the 
fluid left behind the body, and thus a progressive dissipation of tlie lost energy is realized. 

We thus see that the whole question of tlie problem of air resistance consists ui finding 
out the conditions which determine the kind of flow around a body, and we see now how far 
the first attempts to calculate the fluid resistance were from reality. 'Phe.v can only be con- 
sidered as attempts to draw the conclusions from certain assumptions, and it is only with time 
that the idea of the conception of a real fluid, which was always problematic, has slowlj^ been 
reached; and we find ourselves now only at the beginning of tlie development of this great 
question. 

This pampIJet must be considered only as an introduction to tlie question of the law of 
air resistance of aerofoils, which will give a general review of the present main knowledge 
of that question. But a special attempt will be made to show the insufficiency of man\- con- 
ceptions often admitted, and to indicate the ways in which, it seems to me, future investigations 
must be undertaken. At the end I have added some notes which I think will be of interest 
for those who would like to have more complete references concerning the questions discussed. 

Among the questions contained in this pamplilet the following are taken from the "author's 
lectures, given since 1912, at the Polytechnical Institute of Petrograd: The scheme of the 
phenomenon of fluid resistance; calculation of the apparent angle of deflection of the stream 
behind an aerofoil; the establishment of the fundamental wave created by the motion of an 
aerofoil and the determination of its characteristic elements; determination of the part of the 
drag due to tip vorti(;es and its dependence upon aspect ratio; connection between tip and 
edge vortices and the relation of the last to the drag and the lift of the aerofoil; generalization 
of Bernouillis theorem; exact demonstration and generalization of Kutta's theorem; the 
equation of metacentric curves in Flucker's coordinates. 

The author takes pleasure in thanking Dr. J. S. Ames for his kind assistance given by 
reading the manuscript of this Report and correction of its style. 

George de Bothezat. 

Washington, D. C, September, 1918. 



I 



Chapter 1. 



PRELIMINARY CONSIDERATIONS.' 

When, under earth conditions, a certain body is moving, its motion necessarily takes place 
in a fluid, more generally in air- or water. 

When the velocity of the body is relatively low, or the fluid is of low density and low 
viscosity, the action of the fluid medium on the motion of the solid body is not very marked. 
In those cases wo can, without appreciable error, abstract ourselves from the influence of the 
medium and consider the motion of the body as taking place in a vacuum. 

When the velocity of the body reaches a certain value in a viscous fluid of fuiite density, 
the action of the medium on the motion of the body becomes of prune importance. In that 
case, to be able to study the motion of a solid body, we must, in addition to the forces which 
act on the body and among which we necessarily have the Ai-chimedes lift, add a system of 
forces which express the action of the fluid on the different elements of the surface of the body 
in motion. This system of superficial forces, which is distributed over aU the surface of the body, 
is generally called fluid resistance. 

For any body having any general motion in a fluid, the determination of the fluid resistance 
is so complex a problem that its general solution can actually not be found either experimentally 
or theoretically. Only some very simple cases of uniform 
and rectilinear motion of bodies have been, until now, sub- 
mitted to a more or less complete investigation. 

We imagine a solid, which is brought into motion with /• ^ 

a rectilinear and uniform velocity of translation, in a fluid 

medium, which is immobile with respect to the earth, which f 

has uniform and constant temperature, and which has such 
dimensions that the disturbances caused by the motion of the ^' 
solid do not reach the boundary surface of the fluid. In that condition, at a time which is 
generally somewhat after the body has reached this constant velocity, certain steady conditions 
are established. The solid is, so to say, accompanied in its motion through the fluid by a 
certain state of disturbance of the fluid around it. There was a time when it was thought that 
this disturbance has, relative to the body, an invariable configuration; but we now know that, 
generally, this disturbance is invariable relative to the body only before the same, and that 
behind it we often have a state of periodical disturbance. The result of this disturbance is a 
system of steady or periodical forces acting on the whole surface of the solid. This system of 
forces, which constitutes the fluid resistance, can always be brought to a resultant wrench, 
whose comj^ononts will be designated by R, for the resultant force of the wrench, and by C, 
for the resultant torque of the \vTench. (See fig. 1.) If these above-mentioned forces are 
periodical, we will understand by B and G the mean values of the resultant force and the 
resultant torque of the wrench. 

It is quite possible that, for the same body brought into motion with the same velocity, 
the system of forces of air resistance may be difl^ercnt, depending on the manner in which the 
body is brought to its state of motion, but it seems that in most general cases, the viscosity 
tends, so to say, to make uniform all the possible types of disturbances around the body, so 
that, generally, in a free fluid, the same distiu-bances are always established aroimd the body when 
it reaches the same velocity in the same fluid. In that sense we can say: 

For a solid body, moving in a fluid medium with a constant velocity, there corresponds a 
determinate fluid resistance. 

I The main contents of this chapter are taken from the first chapter of the Author's ' ' Etude de la Stability de 1' Aeroplane " Paris 1911, 
a chapter which was written at that time under the influence of the lectures of Paol Painleve. 
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We therefore see that, when the body has reached a steady state of motion, the resultant 
force R and the resultant torque C of the wTcnch of the fluid resistance are independent of the 
time, and are functions only of the magnitude of the velocity and its orientation toward the 
body. 

Under the "law of fluid resistance for unifoim translation of a solid," we wiU understand the 
formulas which give, in position and magnitude, the resultant force R and the resultant torque C 
of the wrench of fluid resistance, as functions of the characteristics of the form and the dimensions 
of the body under consideration, and the magnitude of the velocity of the body relative to the fluid 
and its orientation toward the body. 

It must bo remarked that the components R and C of the resultant wrench do not replace 
fully the fluid resistance, but are equivalent to the system of forces of fluid resistance in only 
one single sense, namely, that modilication which the fluid resistance introduces in the motion 
of the body will be the same when we replace the system of forces of fluid resistance by the re- 
sultant wrench. In all other relations R and (7 are not equivalent to the fluid resistance; for 
example, the stresses which arc produced in the body by the system of fluid resistance are entirely 
modified when we substitute R and C for the fluid resistance. The resultant wrench of a system 
of forces is only an analytical transformation, the possibility of which is established b}- the 

theorems of mechanics and which allows us to reduce a 
given system of forces to its simplest expression. The 
resultant wTench is determined onh' with one degree of 
freedom, its position on its line of du-ection being entirely 
arbitrary. There is no interest in seeking for an exact 
position of the wrench on its line of direction. These 
data would not give us any complementary indication 
upon the motion of a solid. The motion of a solid is 
absolutely determined when the resulting wrench of the 
system of acting forces is given in magnitude, direction, 
and sense. The position of the wrench on its direction 
does not enter into the question of motion. 

It is easy to see that the system of forces of fluid 
resistance can never be reduced to a single resultant 
(orquo, because, if that were possible, the body once brought to that velocity at which this 
could take place would be able to move of itself infinitely forward without any expense of 
power, because it would only be necessary to equilibrate by an acting torque the torque of 
fluid resistance; and this is in full contradiction to all we know about fluid resistance. 

For the same reason the projection of the resultant force R on the direction of the velocity 
must always have a sense inveree to that of the velocity, because if it were not so, the body 
once brought to that state of motion at which that could happen would be able, for example, 
to pull somothmg infhiitely — the torque of fluid resistance being equilibrated by an acting 
torque — and so do work of itself, which would be in contradiction with the principle of energy. 

When the body under consideration has a plane of symmetry parallel to its velocity, the resistance 
of the fluid is reduced to a unique resultant force R lying in the plane of symmetry of the solid and 
whose projection on the direction of the velocity has always the inverse sense of the velocity. 

This proposition can be easily justified. The system of forces of resistance will then be a 
symmetrical system (see fig. 2) and can alwaj's be reduced to a system of forces lying in the plane 
of symmetry; but the latter system of forces can always be reduced to a resultant force or to a 
resultant torque. As we have seen, however, the reduction to a torque being impossible, the 
system of forces will reduce itself to a single resultant force, the projection of which on the 
velocity must have a sense inverse to the velocity for reasons already indicated. 

AH the foregouig does not exclude the possibility of the body's taking a rotary motion 
as a result of the translatory motion in a fluid. 
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We will have especial!}'- to study the law of air resistance of aerofoils. The model aero- 
foil generally has a perimeter of approximately rectangular shape and its cross section is l)\iilt 
up by a system of arcs. h\ figure 3 is represented in plan and in cross section through its plane of 
sjTumetry, an aerofoil of tlie type mentioned, which is movmg in air with a constant velocity 
V pai'allel to its plane of symmetry. Tlie force of air resistance! admits ncccssaril}' a resultant 
force R, whose projection on the direction of the velocity is in the inverse sense of the velocity. 

I once more emphasize the fact tliat no point of the resultant force lias to he distin- 
guished from any otlier point and that the resultant force is fully specified when we know, 
first, its nvagnitudo; second, its line of action, which is defined by direction and position toward 
the aerofoil under consideration. The position of the resultant force on its line of action is 
ad libiium; that is to say, no special point can be, from a meclianieal standpoint, exclusively 
distinguished as pouit of application or center of pressure 
of the resultant force R. 

Very often the couception of point of application of a force is much misunder- 
stood. Let us consider, for example, a nail fixed to a solid body, which is pulled by 
a rope fastened to that nail. In such a case we can certainly speak of the iwint of 
application of a force to our body, which point of application is the xur/ace of contact 
of the nail with the body. When the nail is small and the body large, we can 
abstract ourselves in a first approximation from the size of the surface of nail contact 
and consider that surface as a pphil, and in tHat sense speak of the jioinl of applica- 
tion of the force to the body. But as soon as we begin to speak about the etfuili- 
brium or tlie motion of our body, considering it as a .solid body (by which we meaii 
that we are neglecting the deformation of the body) under the action of the acting 
force, the conception of (X)int of application loses every mechanical sense, as follows 
directly from our statement of the question, because, considering only the equilibrium 
or the motion of the bod}', we abstract ourselves from its other physical properties to 
which belongs also its elasticity. But the consideration of our body as invariable 
brings with it at once that the action of a force upon a body in sense of motion or 
equilibrium is indefiendent of the position of the force on its line of action, so that 
from the standpoint of mechanics of rigid bodies only the magnitude of the forces 
and theirlines of action have to be considered, the position of the forces on their lines 
of action being anyone, and we do not need to consider any point of application. 
Nothing astonishing must be foimd in that last fact. We must only remember the 
whole statement of the problem of motion of the rigid body. It must not be thought 
either that the consideration of such abstract concepts as the rigid body is something 
exclusive. On the contrary, one of the most important scientific methods consists 
in the sorting of the different sides of the questions studied by isolating by abstraction 
a physical property of a body from its other properties; and this general scheme of 
the evolution of our knowledge must never be forgotten, for doing so is the cause of 
great misunderstanding, as happens with the center of pressure in aviation. 

In the beginning of the development of aviation, and by some writers imtll now, 
it was considered as eviden t that the center of pressure, being the point of application 
of the force of air resistance, the airplane had to be considered as suspended at that point when flying. This conception brought at once the false 
conclusion which is in full discordance with experience, that the lowering of the center of gravity would increase the stability, and this false con- 
ception was only the result of the consideration of the notion of point of application, which, exactly speaking, has nothing to do with the motion of 
rigid bodies. By the aid of the theorems of mechanics we can easily Tmd at which point we can consider the airplane as suspended when in flight, 
so that our conclusions fully coincide with reality. This point is the center of m8,ss, because the theorem of moments of momentum is applicable 
to the center of mass indciiendently of its state of motion, so that the oscillation of a rigid body aroimd its center of mass is the same as if the 
center of mass was immovable. From this right conception wo see in full agreement with experience that the weight can have no influence on the 
stability of the airplane, which can be secured only by the forces of air resistance. I have stopped on the last question a little more than I ought 
to do, but the conception of center of pressure is generally so misunderstood in aviation that I thought that this explanation would not be 
unavailing. 

Until now we have admitted that the fluid was immobile and the body moving in the 
fluid, but we could also consider the body as immobile and the fluid ruimmg by the body in 
a uniform stream. If in both cases the realtive velocity of the stream toward the body is the 
same, the flow around the body can be the same in both cases, if the necessary precautions 
are taken for that purpose. But it can very easily happen that ui these two cases the flow 
may be different because of differences in the boundary conditions. 77te piinciple of relulivity 
of hydrodynamics consists m admitting that the fluid resistance depends only upon the relative 
velocity, of the fluid to the body. It is clear that this principle can be admitted only when 
the flow around the body in both cases is the same; and mider the latter conditions the prin- 
ciple of relativity is fully verified and is the (conclusion of the general law of dynamics. 
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Chapter II. 




THE EMPIRICAL LAWS OF AIR RESISTANCE OF AEROFOILS. 

Let us consider an aerofoil represented in cross section by figure 4, moving in air with a 
velocity c(mstant in magnitude and direction. The air resistance B of the aerofoil is fully 
specified by — 

I. Its magnitude; 

II. Its position and orientation toward the aerofoil. 

To find the position and orientation of the air resistance R as well as tlie orientation of 
the relative wind velocity V toward the aerofoil let us take as reference line an arbitrary line 

LL invariably connected with the aerofoil cross section (see 
fig. 4). Wewill designate by a, and call it "angle of attack," 
the acute angle which the velocity V makes with the line 
LL ; by j3 the angle which the air resistance R makes with 
the normal to that same line, and by C the point where 
the line of action of the air resistance R cuts the line LL, 
which point will be called center of pressure. Tlie orienta- 
tion of the aerofoil relative to the velocity F is fully speci- 
fied by the angle a. 

It has been shown by numerous experiments that the 
resultant air resistance encountered by an aerofoil, for certain intervals of the velocity varia- 
tion, follows the following empirical law: 

I. In magnitude the air resistance R of the aerofoil — 

(1) Is proportional to the area A of the aerofoil; 

(2) Is proportional to the square of the velocity Fof the aerofoil relative to the air; 

(3) Is a function of the orientation of the aerofoil toward the relative velocity V; 

(4) Is proportional to the air mass density 5. 

II. In position and direction the air resistance of an aerofoil is independent of the mag- 
nitude of the velocity V and depends only upon the orientation of the aerofoil toward the 
relative velocity. 

The forego mg empuical law of air resistance of aerofoils can be stated in tlie following 
formula : 

R = k5AV-f{a) 

in which i isa coefficient of proportionality and / (a) a function of the angle of attack which 
is characteristic for the type of aerofoil considered. The last formula can also be written: 

R=ICAV'f{a) 

or 

where 

the coefficients Ka and fc<, being certain functions of the angle of attack a only. 

It is customarj' in aerodynamics to consider the resultant air resistance R decomposed 
into two components, the drag Rx, along the relative velocity; .•. 

Rx = R sm + K^A P sm (P + a) 
and the lift Ry, along the normal to the relative velocity. 

R„ = R cos {fi + a) = KcA cos (;8 4- a) 
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when the air resistance i2 is in direction independent of the magnitude of the velocity V, wo 
can write 

where 

K^ = kJ=K„ sin (/3+a) 
X,, = ky5 = Ka cos (/3 + a) 
the coefficients 1cx and being functions of the angle of attack a only. 

Let us now examine the exact meaning of the foregoing empirical laws of air resistance 
and the restrictions to which these laws are sul)mitted. 

We shall begin by an exact statement of the definition of all the quantities which occur 
in the foregoing laws. 

The angle of attack, — Let us designate by a, the angle of attack measured from one ref- 
erence line Z,,Z, and by Oj, the angle of attack measured from a second reference line LJi,^, con- 
nected invariably both with the same aerofoil (see fig. 5) ; and let e be the angle between these 
two lines. It is easy to see that we have 

Oj = a, + € 

If the direction of V varies, we shall have 

+ Attj = «! + Aa, + 1 

or 

= Aa, 

because < is a constant angle. 

We therefore see that the variation of the angle of attack is the same for the same variation 
of the velocity orientation, independently of the reference line from which the angle of attack 




FlO. 5. Fig. 6. 



is measxu-ed. It is probably for the last reason that in the beginning of the development of 
aviation it was thought that the reference lino used to fix the angle of attack can be chosen 
arbitrarih', and the chord of the aerofoil was generally adopted as such reference line. There 
would be nothing to say against such a convention if we had to do only with aerofoils 
with cross-sections of the same type, but all the difficulties begin when we wish to compare 
aerofoils with cross-sections of different profiles. It is in the conception of chord that the whole 
misunderstanding lies. In geometry the word chord is defined as a straight line joining two 
points of a ciurve, but what is the chord of an area like the section of an aerofoil % Nobody knows 
exactly, but, what is still worse, is that it is impossible to establish such a definition. When 
the cross-section profile of the aerofoil is formed by two curves which cut one another, we 
instinctively take as chord the common chord of the curves which limit the profile considered 
(see fig. 6a) ; but for profiles such as represented on fig. 6b two such chords can already be drawn. 
We are still more perplexed for the choice of the chord in the case such as shown in fig. 6c, in 
which any line dra^vn through that profile could with equal success be considered as chord. 
From these simple example.^ we see that the celebrated chord is nothing else than a reference line 
which is chosen arbitrarily. In such conditions when could we say that the profiles a, h, and c 
of the fig. 6 have the same angle of attack.' When we have to do with a flat plate the defini- 
tion of the angle of attack presents no difficulty. It is evidently the angle between the relative 

> This last question is of Grst importance (or aviation practice. For example, how can we judge (or two airplanes having their wings o( different 
cross-section, that they are flying under the same angle o( attackT 
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velocity and the plate itself (see fig. 7). But what is the aerodynamical characteristic of the 
direction of the flat plate? It is nothing else than the direction for which the lift of the plate 
is zero. When the wind blows along the plate, the whole air resistance is reduced to drag 
and we have no lift. It is consecjuently from the direction of zero lift that we measure the 
angle of attack of a flat plate. Thus the direction of zero lift forms our reference line in that case. 

Many years ago Paul Painlcve indicated that, if we wish to obtain a rational basis for the 
establishment of the definitions of all the conceptions which we use in connection with the air 
resistance law of aerofoils, we must simply draw a parallel between the aerofoil and the flat plate 
considered as a conceptional standard. 

Adopting this standpoint, we shall adopt as reference line of each aerofoil the direction for 
which its lift is zero and we shall call that line the zero lift line or, more simple, the zero line. 
The plane normal to the symmetry plane of the aerofoil and containing the zero line will be 
called the zero plane. The zero plane and the zero line are experimentally fully determined for 
each aerofoil. 

Let us consider an aerofoil (see fig. 8) on which the wind blows successively in the directions 
^i) ^'21 ^3) and let Ri, R2, R3, Rt be the air resistance corresponding to those directions. 
We reach the zero line when the resultant air resistance is in the wind direction, as is the case 
for R^. The zero line has to be determined experimentally not only in direction but also in 
exact position relatively to the aerofoil. 

It is easy to see that for each type of aerofoil we generally have four zero lines as shown on 
figure 9. We shall adopt as standard zero line the one which corresponds to zero lift when the 
wind is blowing on the entering edge. It is the zero line which corresponds to T', and Ri in 
figure 9. The angle of attack measured from the standard zero line wiU be designated by i and 
called absolute angle of attack or absolute incidence (see fig. 10); distinguishing this angle from 
the relative angle of attack a measured from any other reference line. The standard reference 
line and absolute angle of attack as above defined are important aerodynamic characteristics 
of the aerofoil.' 

The partisans of chord have reproached the definition of the standard reference line with 
the fact that it is difficult experimentally to measure the incidence from that line. But it is 
quite another question when we have to determine in experimentation the orientation of aerofoils 
in the wind current. In that case we certainly must choose as reference line that line from which 
the measurements are most easily made and such a lino could be called the experimental reference 
line. The question of experimental reference line is a question of the technic of experimentation. 
In one experimental method, one line is more convenient; in another method, another line is 
more convenient. But when stating the results of our experimentation, we must always give 
them in absolute angle of attack, because only in this case will comparison be possible. 

Finally I must also mention the following fact: It can happen that for a certain aerofoil 
cross-section the lift may be zero for any direction of the relative wind within a certain angle, 
as showTi in figure 11. In that case one of the extreme zero lines, V^R^ or V^R^ of the above 
figure, ought to be taken as reference line. In such a case, the lift curve plotted, for example, 
as function of the incidence would have the shape shown in figure 12. 

The aerofoil area. —The area of an aerofoil also needs a special definition. According 
to our standpoint of a parallel drawn between aerofoil and flat plate, we shall adopt as "aerofoil 
area" the area of the projection of the aerofoil on its zero plane. (See fig. 10.) Only with such 
a definition will be avoided all the difficulties and indeterminations, as will be easy to see from 
the detailed discussion which has been made for the angle of attack. 

The center of pressure. — To avoid difficulties, we must also adopt as center of pressure 
the point of intersection of the zero line with the resultant force of air resistance R. 

' For example, an absolute incidence of fivo degrees moans that a five degree decrease of tlie angle of attaclc bring us to icro lift. For the aero- 
foils actually used in aviation practice the standard zero line is generally disposed above the aerofoil, which means that when the absolute 
incidence is equal to zero the air resistance gives rise to a moment relative to the entering edge. 
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I shall show in a few words what, for example, happens if we take for the center of pressure 
the intersection of the so-called chord and the air resistance li. In figure 13 is represented 
an aerofoil and the system of the resultant forces of air resistance Bat ^21 -^4; ^^i' 
different angles of attack increasing in an arithmetical i)rogrcssion. If we follow the dis- 
placement of the center of pressure C referred to the chord, we easily see that there is one posi- 
tion of the air resistance when it is parallel to the chord, and the center of pressure goes to 
infinity. So that the curve of the center of pressure taken on the chord has alwaj's for small 
angles of attack an asymptote. The general shape of the cui-ve of cent«r of pressure in that 
case is shown in figure 14a, where is plotted the distance of the center of pressure 0' from the 
entering edge A' as function of the relative angle of attack a'. If we take the center of pressure 
on the zero line, then the curve of center of pressure will not have any point at infinity and 
the curve of center of pressure will have the shape represented in figure 14b', where is plotted 
the distance AC (see fig. 13) as function of tlic absolute incidence i. The passage of the center 
of pressure to infinity when taken on the chord is only a consequence of a bad definition, because 
it must be remembered that in the definition of the center of pressure we must be guided only 
by convenience. 

To illustrate fully the meaning of the conception of the center of pressure, I shall draw a 
parallel between the notions of center of pressure, center of mass, and metacenier. 

Let us first consider two parallel forces F, and F, of constant magnitude applied at two 
points 1 and 2. (See fig. 15.) As well known, the resultant of these two forces will be 
parallel to them and will divide the distance 1, 2 in inverse ratio to the forces and F,. If 
we consider now the two forces F^ and F^ turning around their points of application but main- 
taining their magnitude and remaining parallel, the resultant force Iti^ will also turn around 
a definite point. If we consider now a system of three parallel forces F„ F^, F^ constant in 
magnitude and turning around then- points of application, it will be easUy seen that the 
resultant force of the thi'ee forces will also turn around a definite point, because i?^, is the 
resultant of and F^, and so on, independently of the number of forces. The point through 
which the resultant force of a system of constant parallel forces turning around their points 
of application always passes, is called the center of the parallel forces. The center of mass is 
a particular case of center of parallel forces when the forces considered are the weights of the 
different elements of a body. 

Let us now consider generall}^ anj- system of forces applied at any points. If we con- 
sider the continuous variation of these forces, their resultant force will also vary continuously in 
magnitude, position, and direction, and will describe in space a certain surface which is called 
the mttacentric surface. When all the forces considered lie in the same i)lane, the resultant 
force also lies in the same plane and the metacentric surface is reduced to a metacentric curve, 
which is the envelope of the successive positions of the resultant force. The point at which 
the resultant force touches the metacentric curve is called the metacenter. (See fig. 16.) When 
the forces considered are parallel and constant in magnitude, the metacentric curve reduces 
to a point. We therefore see that we can consider the center of mass as a particular case of 
metacentric curve reduced to a point. 

If we consider the system of forces of air resistance of an aerofoil, these forces admit a 
metacentric curve and it will be easy to see that this metacentric curve has alwaj's a cusp 
point admitting the zero line as tangent at that point. In figure 17 is represented the general 
shape of the metacentric curve of an aerofoil. For comparison, in figure 18 is represented 
the metacentric curve of a fiat plate.' 

The important fact is that the center of pressure is neither a center of parallel forces nor 
a metacenter, but simply a point arbitrarily chosen to fLx the position of the resultant force 

• The scale used on fig. 14a to plot the distance A' C is smaller than the scale used on fig. lib to plot AC. 

' In the theory of the airplane the metacentric curves do not have the same importance as in ship theory. In the last theory the metacentric 
curves allow a direi:t evaluation of the restoring moments, on account of that fact that, to a first approximation, the lifting force of a ship is con- 
stant when the ship undergoes oscillations. It is not the case for the airplane, where the lifting forces are variable in magnitude when the airplane 
is oscillating, so that the metacentric curve alone docs not determine the restoring moment. That is why for aerofoils the metacentric curve must 
be considered only as giving the general picture of variation of the resultant force of air rcsitance in position and direction. 
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of air resistance. That is why we must choose it in the way most convenient for our purpose. 
For this last reason we shall take the center of pressure on the zero line. In that case we shall 
have no point of the center of pressure curve in infinity, and the cusp point of the meta- 
centric curve will be the limit position of the center of pressure on our zero line. 

The velocity. — The velocity F which enters in the formulas of air resistance of aerofoils 
disposed in a uniform, fluid current has to be taken in front of the aerofoil and at such a dis- 
tance from it that tlie disturbances in the medium caused by the presence of the aerofoil do 
not reach it. This distance generally lies in front of the aerofoil between one and two 
times its breadth. 

The fictitious equivalent plane. — As a summary of all the foregoing discussion we are 
brought to the following conception: 

Let us consider the zero plane of an aerofoil and project on this plane the aerofoil and 
take this area as the jiditiom equivalent plane, or, shorter, as equivalent flane, of our aerofoil; 
that is to say, attribute all the properties of our aerofoil to that plane and refer all the quan- 
tities which we use to describe the law of air resistance of aerofoils to that fictitious equivalent 
plane. We shall thus take as area of the aerofoil the area of the equivalent plane. (See fig. 
19.) We shall measure the angle of attack from that equivalent plane and this will be our 
absolute incidence i. We shall take the center of pressure on that equivalent piano and fix 
the direction of the force of air resistance R by the angle /3 of its inclination to the normal to 
that equivalent plane. Under such conditions, all the formulas of pages 2 and 3 have to be 
referred to the equivalent |)lanp: and in that case we shall write: 

I. The magnitude of the air resistance R of aerofoils (see fig. 19) 

R = ^5^ V'f{i) = liA V'fd) 
R=KiAV^ = 1cMV^ 

where 

K=k5; I{i=Kf(i)=k^ 

The drag 

R^^Rsin. (/3+i) = /4^P 

The lift 

Ry=R cos W+i) = KyAV' 

where 

Ky = Jcyd = Ei cos (jS + i) 

the coefficients K^, K„ and k^, l'„ being functions only of the absolute angle of attack i. 

Some general data on aerofoils.— To the foregoing formulas I will add the following 
remarks: 

For the orientation under which the aerofoil is practicaDy used, the lift of the aerofoil 
is generally equal to zero only when the wind is blowing on the back of the aerofoil, and the 
equivalent plane is disposed somewhat above the aerofoil. The position and orientation of 
the equivalent plane can, in general, also depend from the value of the speed V, so that to 
different speed intervals can correspond, for the same aerofoil, different equivalent planes. 

Starting from zero absolute incidence, the air resistance R rises very quickly out of the 
zero plane, so that for angles of attack around 5°, the air resistance makes small angles with 
the normal to the zero plane. 

For the aerofoils actually used in aviation for small angles of attack, the ratio of dragf 
to lift can reach 1/20. 

For actual aerofoils, considering the incidence increasing from zero, the center of pressure 
first approaches the leading edge (sec fig. 17) — that is, travels in a sense inverse to that for the 
case of a flat plate (see fig. 18) — and only afterwards, for greater values of the angles of attack 
(generally larger than 10°) the center of pressure begins to travel away from the leading edge. 

Tlie coefiicients Kx, Ky, and k^, k„, for equal values of the angle of attack, have the same 
values only for aerofoils having similar cross-section and similar perimeters; and still in that 



= kMV' 
^kybAV 
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case their values also depend upon the magnitude of the aerofoil area A and the magnitude 
of the velocity Y. These coefficients do not vary much when the magnitude of the aerofoil 
area changes, and ])y the variation of the velocitj^ the coefficients K^- and /iv are principally 
affected — they diminish when the velocity increases — the coefficients Ky and k„ do not seem to 
depend much upon velocity for a value of the last above a certain value. 

For small angles of attack — up to around 10° — for most actual aerofoils, the coefficients 
follows a linear law and the coefficient a parabolic law. So that for such angles we can 
write 

]cx = Jc{ai' + hi + c): 

so that in such a case the lift and drag of the aerofoil have for expressions 

Ry = UAV^i= KAVH 

R^^UAViai' + hi + c) = laVHai^ + bi + c) 

The value of the coefficient K depends upon the aspect ratio L/b, that is, the ratio of its span L 

to its breadth b. For values of this ratio equal to about 
5 or more, the coefficient K for most actual aerofoils, for 
usual atmospheric conditions, has a value near to 1/200, 
the units used being the meter, the kilogram, and the 
second. For smaller values of the aspect ratio, the 
value of K diminishes. 

Different characteristic curves used to plot 
the results of measurements of the air resistance 
of aerofoils. — To plot the results of measurements of 
air resistance of aerofoils different systems of curves are 
used. From any system of characteristic curves giving 
a full specification of the laws of air resistance of aero- 
foils, we can deduce any other one. 

First method} — The most direct way of representing 
the air resistance of an aerofoil is to plot the curves of 
the coefficients Ki or Z\ as function of the angle of 
attack i, and the curve of the angle /3 as function of the 
angle of attack i. The Ki or /(•,■ curve gives a direct 
evaluation of the magnitude of the force of air resistance; 
and the /3 curve gives the laws of variation of the incli- 
nation of the air resistance to the normal to the zero fine. The general shape of the Ki and /3 
curves are represented in figure 20. 

Second method. — Another method very widespread in the practice of modern aerody- 
namical laboratories is to plot the lift curve Ky and the drag curve Kx as functions of the 
angle of attack i. To these curves the drag-hft 7ix/A"„ curve is generally added. It is much 
more convenient to use the drag-fift curve than the lift-drag curve, as it is made sometimes, 
because many fundamental properties of the airplane are directly connected with the drag- 
lift curve.^ The general shapes of the drag curve, the lift curve, and the drag-lift curve are 
represented in figure 21. 

Third method. — Probably one of the oldest methods used to represent the laws of air 
resistance of aerofoils consists in plotting the lift coefficient as function of the drag coefficient. 
This method was used by Lilienthal. When using this method the angle of attack is marked 
on the curve. (See fig. 22.) This method presents the advantage that we can also read on 
the curve Zy= F{Kx) the variation of the drag-lift ratio. It is easy to see that the tangent 

' This method of plotting is used by the author for propeller calculations. 

• Sec, for example, G. do Bothczat, "Etude dc la Stabihlic de rAcroplane," pp. 58-68, Paris, 1911. 
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of the angle 7 which a hne joining the origin with a point of the curve /r„= FiK^) makes with 
the Ky axis is equal to 

so that we can directly plot a scale for Kx/ K,, on a parallel to the Ki axis. Each straight 
line joining the origin with a point of the K,, = FiK^) curv-e cuts off on that scale the value 
of K-c/K,,. The tangent drawn from the origin to the curve Ky= F(,Kx) gives the minimum 
of the value of Kxl Kv 

In his last research on aerofoils, Eiffel uses this method and for convenience plots the 
at the scale ten times bigger than the Ky. 
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Fig. m. 



To specify fully an aerofoil by each system of the foregoing curves, there must be added 
the curve of the center of pressure and the zero line in exact position and direction. It is 
also good to draw the metacentric curve which gives a full picture of the positional and direc- 
tional variation of the forces of air resistance. 

I must also add that it is necessary that the data on aerofoils be at least determined for 
an inten^al of —90° to 90° of absolute incidence. Tliis is on account of the fact that we must 
not limit ourselves to the actual necessities, but Jiiust also give data which future research 
and discovery may need.' 

' As example, I can indicate tlie (ollowinc (act: Only Iwcause the aerofoil data were not enough extended, we can not actually calculate the 
thrust of a propulsive screw at a lixed point whenever the necessary methods iire already at our disposal. A propeller at a fixed point works at 
very large angles of attack, 20°, 30°, 40°, and in some cases still greater. 
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Chapter III. 



THE FLOW AROUND THE AEROFOIL. 




Fig. 23. 



In the begimiuig of the development of aviation, the maui knowledge arose at first only 
to the general quantitative relations about air resistance of aerofoils, and it is only veiy slowly 
that light has been throwni on the (low phenomena. 

For an aerofoil moving with a uniform and rectiluiear velocity in air the following flow 
phenomena have been progressiveh^ discovered experimentally. 

I. On the ui)per surface of the aerofoil we have a decrease of pressure and on the lower 
surface of the aerofoil we have an increase of pressure. Tlie depression created 07i the upper 
surface is, for small angles of attack, always larger than the increase of pressure on the lower 

surface, so that the lijt of the aerofoil is due 
more to a sxtction exerted on the tipper side than 
to the pressure exerted on the lower side. 

II. The stream in the wake behind the 
aerofoil appears to be deflected do^vniward. 

III. From the tips of the aerofoil vortices 
run off which we wUl call the tip voiiices (see 
Fig. 23). The rotation of the fluid hi these 
tip vortices has the sense fi'om the outside 
space into tlic inside space between the vor- 
tices, if we look from above.' 

IV. In the space between the tip vortices 
two kinds of flow can take place. For very 
small angles the flow is contumous; that is to 
say, we have no sensible turbulent motion. 

But when the angle of attack increases beyond a certain value of the last, there appear on both 
edges of the aerofoil vortices, parallel to these edges, which we will call the edge vortices. The 
greater the velocity of the flow running on the aerofoil, the smaller is the angle of at tack for which 
edge vortices appear. These edge vortices are not stationary with reference to the aerofoil. 
They grow up on the edges of the aerofoil and, when they have reached a certain intensity, they 
nm off ui the general direction of the stream behind the aerofoil, so that these edge vortices 
have a certain velocity with reference to the aerofoil. 

The edge vortices which grow on the upper and lower edge rotate from the space outside 
the two vortices into the space inside, when one looks from above, so that behind the aerofoil 
there appears a system of vortices in quincimx ^ rotathig in inverse senses (see fig. 24). The 
ends of these vortices go over into flio tip vortices. As the edge vortices are rotatbig in inverse 
sense the mean value of the intensity of the tip vortices is not modified by the edge vortices. 
We now see that the general picture of the flow behind an aerofoil looks like a vortex ladder 
running off the aerofoil.' 

So far as I know, exact measurements of the depression on the upper side and the pressure 
on the lower side of an aerofoil were firet made by G. Eiffel. The apparent stream deflection 
behin<l the aerofoil seems to have been observed by many investigators. Tlic ncceasitj' of the 
existence of tip vortices seems to have been fu-st indicated by Lanchester. The vortices in 

' Observations of tip vortices have been made by J. R. I'annolland N. R. Campbell. "Tiie Flow of .\ir Around a Wing Tip," Report No. 197, 
Marcli, 1916. Advisory Committee for .Aeronautics. 

' I call7«i'>c'"!i voHa 5»»/cn). a system of two parallel rows of equidistant rectilinear and parallel vortices rotating in inverse senses in each row, 
and in such an arrangement that the vortices of one row are disposed towards the middle of the distance between the vortices of tlie other row. 

' Good pictures of edge vortices running off an aerofoil can be found in the " Teclniical Report of the .Vdvisory I'omniittee for .\eronautics," for 
the years 1912-1913. " t'hotograpbic Investigation of the Flow Around a Model Aerofoil " by E. K. Rclf. p. 13.1. 
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quincimx were first noted by Karman for the particular case of the orthogonal motion of a flat 
plate and the motion of a cylinder. 

Before attacking the detailed discussion of the above-mentioned flow phenomena, I shall 
first make some general remarks on the flow phenomenon aroimd a solid body moving in a fluid 
with a uniform velocity. 

General scheme of the phenomenon of fluid resistance. — I shall here develop a general 
scheme of the fluid resistance phenomenon, which must be considered as a conceptional limit, 
but which takes account of all the fundamental circumstances of the fluid resistance phenomenon 
in their most simplified form, and thus allows a better view of the relations which hold. 

We shall first consider the case of an infinite cylindrical body having a plane of symmetry 
and moving in an infinite fluid with a constant velocity Fo parallel to that plane of symmetry, 
the sense of Fo being taken as positive sense. (See fig. 25.) 

Let us imagine an observer moving with the bodj^. For such an observer there will appear 
a relative stream runnuig on the body. 



0^ 




Fig. 24. 



When the fluid is considered as perfect this 
relative stream can be assumed as benig a 
potential stream — that is, a stream ad- 
mitting a velocity potential for the velocity 
distribution in it. But for a real fluid, in . 
the case of our problem, according to the 
indications of the experiment, there must 
necessarily be losses inside the fluid and 
thus a certain distribution of vortices in it. 
This last fact is a direct consequence of 
the general equations of motion of a viscous 
fluid, according to which there are no losses 
inside the fluid where there are no vortices.' 
For the general analysis of the fluid resist- 
ance phenomenon we will place ourselves in 
ideal limiting conditions and replace the 
effective relative stream rmming on the 
body by a conventional relative stream, but 
so defined that in relation to the fluid resistance the conventional relative stream will be fully 
equivalent to the effective relative stream. 

We will first assume that in each cross section normal to the plane of symmetry of our body 
the velocity of the conventional relative stream is constant. In such conditions, to take 
account of the change in the distribution of the velocities in the general stream which are pro- 
duced by the presence of the body, we must consider our conventional relative stream as limited 
by surfaces of discontinuity outside which the general stream velocity is uinnodified, but inside 
which the velocity, being constant in each cross section, is different from the outside velocity. 
(See fig. 25.) These surfaces of discontinuity must thus necessarily be constituted by vortex 
sheets. On the other hand, as we must also conceive the fluid as adliering to the surface of the 
body — a fact to which seem to lead Zahm's ^ experiments on skin friction, which have shown 
its independence of the state of the body's sm-face — we must consider the surface of the body 
as covered by a vortex sheet in which gliding of the fluid takes place, the relative velocity at 
the surface of the body being equal to zero. We thus see that in our conception of the flow 
phenomenon around the body the vortices instead of being spread in a certain way inside the 

> According to Lamb, " Treatise on the Mathematical Theory of the Motion of Fluids," the dissipation of onergs' inside a fluid mass is given 
by the general expression 



4 ''J' J" J" ("«'+"y'+"«') (fr i'J <l2 



which is equal to wo for ui— u,=-<j,— 0 

" "Atmospheric Fricliou on Even Surfaces," I'hilosophical Magazine, July, 19tH. 
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fluid have to be conceived as concentrated on the surface of the body and on the boundary of 
the conventional relative stream. (See fig. 25.) 

Let us consider now two cross sections, I and II, of the rehitive stream running on the 
body. Both cross sections are considered immobile relative to the body and are taken normal 
to the stream. Section I, of conventional height ho, is taken before the body at a distance not 
reached by the disturbance created by the body in the fluid. The relative velocity and the 
pressure in that section, uniform in the whole section, are designated by 2'o and — Fo, this last 
velocity being equal in magnitude to the velocity of the body, but having an inverse sense. 
The Section II is taken behind the body; p and — F are the uniform pressure and relative 
velocity in that section. The absolute velocity w of the stream behind the body is equal to 

(1) w=- F-(-F„) = Fo- F 

The velocity w is nothing but the mean velocity of the wake behind the body. 




Fig. 25. 



Let us designate by O the value of the Bemouilli constant corresponding to Section I. 
According to BernouiUi's theorem we must have 

SF' 

(2) Po+'-^=C. 

Let us follow from Section I to section II a stream line in the relative motion of the fluid 
toward the bodj-. Starting from the values po and Fo in Section I, pressure and velocity will 
vary along the stream line. When we pass by the body, the velocity will be increased and 
the pressure decreased. Behind the body the velocity will drop and the pressure increase, 
and when we reach Section II we shall find there a pressure p and a velocity — F connected by 
the relation : 

SF^ 

(3) p+'-~=C-AC, 

where AC is the drop in the Bemouilli constant, which occurs when we go from Section I to 
Section II. In reality this drop is due to the losses taking place at the surface of the body by 
skin friction and inside the fluid by viscosity, which losses, in our limited conception, are assumed 
to be concentrated on the boundaries of our conventional stream. Subtracting equation (4) 
from equation (3), we obtain the relation connecting po and Fo in Section I with p and F in 
Section II. 

(4) Po-?>+|(f„^-F^)=AC. 

Let us designate by n the whole amount of work done by the forces of viscosity inside 
the fluid between the Sections I and II. We can always consider this interior work referred 
to the velocity Fo in Section I and consequently write 

(5) n = FFo, 
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where F is a fictitious force, which we call the dissipative force and which corresponds to the 
interior work done between the Sections I and II and referred to the velocity Vo. By its direct 
meaning, AC is the interior work by unit of volume of the relative flow crossing the two sec- 
tions I and 11 in a second, which volume, counted per unit of length of the body, as a conse- 
quence of continuity, is equal to either of the two expressions 

(6) KV^^hV, 

when one neglects the very small density variations. We thus have 

n FVo F 

or 

(8) F=hJ^C. 

The general picture of the flow around the body begins now to appear more clearly. Outside 
the boundaries SS and S'S' of our conventional relative stream (see fig. 25), as we consider 
the fluid unaffected by the motion of the body, we thus have the uniform pressure po', and in 
the absolute motion the fluid is immobile, so that in Section II from outside to inside we have 
a difi'crcnce of pressures Po — p and a difference of velocities w= Vq— V, maintained by the 
conventional boundaries SS and S'S', which are vortex sheets. 

Let us now apply the theorem of momentum to the fluid mass contained between the body 
considered, the vortex sheets SS and S'S', and the Sections I and II and included between 
two planes normal to the body at a unit distance from one another. Let us designate b}' Rj^ = 
kjbh T o" the drag of the body counted per unit of length, h being a linear dimension of the body. 
This drag constitutes, in our case, the resultant of all the forces acting on the surface of the 
bodj'. It will be easy to see that we have 

(9) - P + UPo -V)= F„H*) 
This last relation in connection with the relations 

(10) (Po-p) + |(TV-T'')=AC=£ 
and 

(11) V^^VTi 

gives thus three equations connecting the pressure p, the velocity V and height A in the section 
II with the corresponding quantities p„, V^, and Ti„ in the section I. From all these quantities 
the only ones to be considered as known in our problem are and r„. 

It is easy to show that when the section II is considered taken at such a distance from the 
body that either p = Pa or F= Fj, we will have very approximately • 

(12) F^h^bh Fo^ 

That is, the dissipative force equals the drag. 

For when F= F„, which brings with it A = A„, we have 

(13) F=\{.p,-p) = Upo-p)='lc^^hV'' 
and when p = Po, then 

(14) F=h,{V^- V'} 



2 

2 ' • 2 



* To And the resultant of the outside pressure on the boundary of the portion of the relative stream considered, it is sulTicient to conceive the 
pressure p» added and subtracted in the section II. We then will have a uniform pressure p, on all the boundary surface, whose resultant is lero, 
and the pressure h (.pt-p) in the section II, which quantity constitutes the resultant pressure. 
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But, as will be seen in the following, when is a small quantity, the difference between h,,, h, 
and Fo, V in case of p = Po is negligible, so that 

(15) F^5h„V,^-5hV' = JcJlV„' 
In the general case 

(16) F<k^dbV„' 

The relations (15) and (16) lead us to make the assumption that the dissipative force can 
be considered as having the form 

(17) F=f5bV,' 

where / is a characteristic coefhcient depending upon the form of the body and the properties 
of the fluid and the position of section II. In the case of p = p^ or V= we have 

But in the general case we will have 

f<K 

The inferior limit for / is determined by the skin friction at the surface of the body.' 

All the foregoing constitutes, so to speak, a limited scheme of the fluid resistance phenom- 
enon, but one which gives a complete picture of the relations occurring. Let us examine the 
connection between our scheme and rejility. 

In the relative flow around a body, the observed velocities in a section such as section I 
are uniform when the section is taken at a suflicient distance from the body, but in a section 
such as II the uniformity of velocity is generally not observed. This last fact does not consti- 
tute an essential difference, because we can always conceive the velocity F as a certain mean 
value of the real velocities. 

Much more essential is the question of the practical possibility of the existence of the 
vortex sheets at the surface of the body and on the boundary of the stream. It has been 
pointed out by many investigators that vortex sheets in viscous fluids must be considered as 
unstable.- Experiments performed on the observation of the flow around bodies, although 
not very numerous, have already given valuable indications.^ For relative flow velocities 
having a sufficient value, the vortex sheet covering the surface of the body always passes over 
into a system of vortices in quincunx. This last fact was first fully understood by Karman, 
who also indicated the reason why we get the quincunx vortex system. Karman's investiga- 
tions of the quincunx vortex system have shown that this system is stable. The edge vor- 
tices above mentioned are nothing else than the vortices i/n quincunx into which the vortex 
sheet covering the surface of the body passes. For low velocities we also have in all prob- 
ability a tendency toward the formation of the vortices in quincunx, but the energy in the 
wake being small, the energy of the beginning vortices is dissipated before their complete for- 
mation. The motion which is established must be a kind of turbulence which distributes 
inside the fluid the vortex sheets covering the surface of the body. The mechanism of this 
distribution is in all probability the foUowing: We either have a direct, irregular, and periodical 
transformation of the surface vortex sheet in quincunx vortices, dissipated before full forma- 
tion, reformation of the surface vortex sheet, and so on; or we have a periodical irregular for- 
mation of the surface of discontinuity established in the Kirchhoff-Lord Raleigh theory. 
These surfaces of discontinuity, which must necessarily be vortex sheets, can appear as inside 
boundaries of the relative motion, only as momentary phenomenon. At such a moment the 
flow appears as represented in figure 26. But these surfaces, being unstable, quickly disappear 
and the vortex intensity concentrated in them is dissipated before the formation of a definite 

> It must be remarked that the calculation of the resultant skin friction at the surface of a body often presents some ambiguity, the exact dis- 
tribution of the velocity at the surface of the body having to be known. 

» Sco H. Poincare, "Theories des Tourbillons," p. 174, § 142; A. B. Basset, "A Treatise on Hydrodynamics," Volume U, p. 309; and C. Schaefer, 
" Einfiihrung in die Theoretische Physik," Band I, p. 896, J 199. 

< For a review of these observations see W. L. Cowley and H. Levy, "Aeronautics in Theory and Experiment," Chap. H. 
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kind of vortex motion. Such irregular, unstable, periodical process, affected by the smallest 
perturbation, spreads the surface vortices in the fluid. 

In wliole probability all the observations made upon the vortex sheets surrounding the 
body apply also to the vortex sheets SS and S'S', constituting the outside boundaries of our 
conventional relative stream. For small relative flow velocities, those vortex sheets are in a 
certain way distributed in the fluid ; but for greater velocities it is possible that they go over into 
the stable system of vortices in quincunx. We thus see that behind a body moving in a fluid 
we shall have, in general, a periodical fluid motion. I shall call primary or fundamental wave the 
fluid motion generated by the formation of quincunx vortices from the vortex sheets limiting the 




FiQ. 28. 



relative stream boundaries from the outside, and secondary wave the fluid motion generated by 
the formation of quincunx vortices from the vortex sheets limiting the inside relative stream 
boundaries. Both fundamental and secondary waves will be considered in more detail in the 
following. Tlie possibility for the existence of the fundamental wave will appear with still 
more evidence from the general examination of the flow around an aerofoil, to which we shall 
now pass. 

The same scheme which we have developed for a symmetrical body can be applied to an 
asymmetrical body like an aerofoil. All that has been said relative a symmetrical body has to 




Fig. 27. 



be directly transferred to the aerofoil. The difference will consist in the fact that, as the aerofoil 
has a lift component due to the fluid, there must be a fluid momentum corresponding to that 
lift. That is to saj', the relative flow beliind the aerofoil must be deflected dowTiward. The 
schematical flow around an aerofoil is represented in figure 27. Let us now imagine for one 
moment the aerofoil immobile and the stream running on it with the velocity Vo- The fluid 
velocity outside the stream boundaries SS and S'S' wUl also bo Vg. In such a condition it is 
easy to see that there wiU be a tendency to straighten the deflected stream by the outside 
stream. If we assume the possibilit}' for the stream between the boundaries SS and S'S' to 
become horizontal after section II, the application of the momentum theorem for the lift, 



\ 
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between two sections such as I and II, will show an increased pressure above and a decreased 
pressure p, below. In such a condition it would be difficult to imagine how the theorem of 
moments of momentum applied to the stream portion between SS and S'S' on one side and I 
and II on the other side coidd be satisfied. Such a flow appears impossible, and it is easy to 
see that, after its downward deflection, the stream, by the difference of pressure (pz-Pi) must 
necessarily be deflected upwards. It thus becomes evident that behind an aerofoil we have 
a wave motion of the relative stream. The instability of the vortex boundary sheets SS and 
S'S' also lead to this conclusion. The wave motion which is to be expected is represented 
schematically in figure 28. It is sufficient to look at this last figure to see at once that the 
wave motion obtained is governed by a system of vortices in quincunx, rotating in one sense 
for the upper row and in an inverse sense for the lower row. We thus see that behind the 
aerofoil we can expect to see the phenomenon of the fundamental wave mentioned in the 
foregoing. 

Tlie phenomenon of the secondary wave can also take the place for the aerofoil. 

We arc thus brought to the conclusion that a simple deflection of the relative stream behind 
the aerofoil is not to be expected. Nevertheless, the preliminary study of the stream deflection 
behind an aerofoil is of interest for the following reason. 




Fig. 2S. 



Let us consider an aerofoil II disposed in the wake of another I. This aerofoil II will 
then be submitted to a periodical stream. Let us assume for simphcity that both the magnitude 
of the velocity V of the flow running on the aerofoil II and the angle of attack % vary according 
to sinusoidal laws, so that as a first approximation we consider 

V= Vo+v sva.-jrt ; i = %o — '] sm^ t 

T being the period, v and j the amplitudes of variation of V and i, t the variable time. The 
difference of sign in the above expression denotes the fact that the velocity V is assumed 
increasing when the angle of attack decreases, and vice versa. We shall also assume that for 
the instantaneous values of V and i the resultant air-resistance R of the aerofoil II can be 
expressed by the formula : 

R=KAVH^\VH 

writing ZA = X. Let us now, under these assumptions, calculate the mean value Rm of R. 
We have 



fi«=4 I Rdt 

)dt 



4jr 



XFo'io / i , XV.\2vi.-jV.) I ''^tdt + 
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As is well known, the integrals of the odd exponents of the Sine are equal to zero, and 
So that we finally get 

For small values of v the term i \vHo can be neglected, and thus 

and we are brought to the following conclusion: It is easy to represent such a periodical flow 
that an aerofoil disposed in it wUl show, for the mean value of the air resistance, an apparent 
decrease of the angle of attack. Tlio aerofoil will thus appear as if placed in a downward 
deflected stream. The study of the apparent stream deflection behind the aerofoil is thus 
justified. 




FIG. 29. 

The apparent stream deflection. — If we make the assumption that all the lift of an 
aerofoil is due to the momentum created by the deflection of the stream downward, the angle of 
deflection of the stream can be easily calculated. This calculation will also give us a mean 
value of the height of the stream disturbed by the presence of an aerofoH in the fluid. 

Let us consider a unit of length of an aerofoil and draw aroimd it a contour a, b, c, d defined 
as follows (see fig. 29) : 

The side ab of the contour is a plane cross section of the stream taken at such a distance 
before the aerofoil that the flow in that section is not disturbed by the presence of the aerofoil. 
It is the cross section I of the relative stream. The two sides ac and hd of our contour are 
taken along two stream fines at such a distance that the local phenomena created by the presence 
of the aerofoil in the fluid do not reach them, and that the pressure on these two stream fines is 
equal to the outside pressure. At the end, the side cd, constituting the section II of the relative 
stream, is a plane cross section taken at such a place that the velocity of the stream has nearly 
taken its original value and the pressure has also nearly reached its original value. These 
last assumptions are only a certain approximation. Let us now calculate the increments of 
the components of the momentum of the fluid running out of this contour for a unit of length 
of the aerofoil, along the velocity and along the normal to the velocity. 

It is easy to see that the increment of the components of the fluid momentum along the 
velocity is equal to 

^6P-^Pcos a = ^y*(l-cos a) 

and that the increment of the components of the fluid momentum along the normal to the 
velocity is equal to 

iWF»sin a 

147498— 20— No. 28 3 
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As these two components of the momentum must be equal to the two components of the 
fluid pressure on our aerofoil, we must have 

(18) 5^ = M6^ = ^^^(1-C08 a) 

(19) Ry = IcySb V = U F'sin a 

where b is the breadth of the aerofoil and S the air mass density. From the last relations it 
follows that 

Arj.6 = A(l — cosa) ; lcyb = h,s\na 



Tc^ 1 — cos « _ . « 
Tcy~ sin a ~ " 2 



but, as we have 



we easily get the value of 



2tg| 

2 



Having the value of sin a we easily get the value of the stream height disturbed by the aerofoil: 



^lCyb_ J ^ + ]^ )_ b iK' + V) 



(21) ^ sin a 2Tc^ ' ' 



For small values of the angle a and negligible values of Tc^, we can wTite 

(22) ^~2t^ 

(23) h^^vL 



The value of these formuhis, the deduction of which is only based on certain assumptions, 
lies in the fact that they indicate from which quantity depends the apparent angle of deflec- 
tion of the stream and the height of the fluid stream disturbed by the presence of the aerofoil. 
Further, these formulas give values for both quantities of the order of magnitude as obtained 
from experiment.' 

The formula (22) is capable of an interesting geometrical interpretation. If we draw a 
plane normal to the resultant air resistance B, we then see that the stream is, so to speak, 
reflected on this plane. (See fig. 30.) 

From these very simple considerations we see that the conceptions developed in turbine 
theoiy, where it has been assumed that the fluid nms off from a turbine wing in the direction 
of the tangent to its trailing edge, are absolutely inadmissable. The apparent direction of the 
stream behind an aerofoil or a turbine wing depends not only upon the direction of the tangent 

> For example, let as take (the units used being m.,kg., and sec.) 

i-6» 

I k,i Ki 1/200-6.^, 

''--T-T J * 

(r,/ll:,-l/16 

we then have 

* 2i. TIT rjm 



'F, 



* 



a quantity whioh Is of the order o( what experience with biplanes indicates to be negligible influence of the mutual Interfetence of the wings. 
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to the traihng edge but upon the form of all the parts of the aerofoil or the wing. The general 
results to which we are brought in this elementary calculation are certainly only of a first ap- 
proximation, but they give a rational description of the general phenomena. 

We shall now proceed to examine the problem of the apparent stream deflection behind 
the aerofoil to a second approximation. 

When the stream meets the aerofoil, as the result of the impact which takes place there 
must be a certain amount of energy dissipated inside the fluid. Let lis designate for the section 
I before the aerofoil, by p„, Vg and ho, the pressure, the flow velocity and the relative stream 
height distiirbed by the aerofoil, and by p, V and 1i, the values of these same quantities for 
the section II behind the aerofoil. We shall apply to the fluid between the sections I and II 
the momentum theorem taking account of the dissipation of energy by two limiting assump- 
tions. The first assumption will consist in considering in the section II the velocity V= V^; 
the second, in admitting p = Po- Under such conditions the dissipative force F will be equal 
to the drag of the aerofoil for an angle of attack equal to zero, as follows from the foregoing. 
For other values of the angle of attack, the coefficient/has to be considered, for a given aerofoil, 
and as a function of the angle of attack. In the case of F= ¥„, as we have F=h{po — p), the 
dissipative force Fc&n be conceived as applied in the section II normally to that section. 




Fio. 30. 



Applying the momentum theorem, in the case of the first assumption, it is easy to see that 
we have: (See fig. 27, and compare with the similar equations (18) and (19) ). 

(24) Rx = h^b F = F - h5 V cos a +f5b V cos a 

(25) Ry = ICySb F = hS F sin a -fSb V- sin a 

which equations express the fact that the drag and lift are equal to the corresponding com- 
ponents of the variation of the fluid momentum to which are added the components of the dis- 
sipative force F. 

Dividing the last equations by 6b F we get 

(26) ix^^/b — h/b cos a+f cos a 
(27) 

from which follows 



lCy = h/b-f sin a 



sin «= 



___ and cos « = i 



and, since sin' a+cos'a= 1, we have 

(A/6-/)'=V+W*-*x)' 
or, removing the parentheses and midtiplying, we finally find 



(28) 

or 

(29) 



2(fe,-/) 
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By aid of the relation (28) we find 

Till- f = V 

and consequently 

Neglecting the squares of fc^ and/ which are of the same order and very small quantities for 
actual aerofoils and small angles of attack, we get 

(32) sin a^a^?^^^pfl 

Ky 

In the case of the second assumption, that is, for p = Po and AC=8/S{Vg^- F*), we have 
(see fig. 27 and compare with the similar equations (18) and (19) ) 

(33) Tc^dhVo' =hoS Vo' -MV'coa a 

(34) Icp Sb Vo' = hS Psin a 
or 

V 

(35) K = K/i- ii/S-fTjCOS a 

' o 

(36) fc„ = A/6^sina 

' 0 

and taking accoimt of the condition of continuity hoVo^liV vre get 

y 

(37) Jiji-Tcx^ hoi 6 T^cos a 

'o •- 

V 

(38) Jcy = holhy^sm a 

' o 

Taking into account the relations (8), (15), and (16), we find 

yi 2F _2fl 

so that the ratio FV Vg^ is seen to be equal to 

V* i 

(39) V?='-'fK 

Squaring the relations (37) and (38), adding them and substituting in the last the foregoing 
value of T^/FoSwefind 

(40) {Klb-hy + Tcy'= Qijmi ~2fhlK); 
from which relation we directly find 

(41) ^olh-l^f^ 

Introducing this last value of hjl) in the equation (38) and substituting in it for F/ Vg its value 
we find 

(^^) ^""^ V(57W)WTit7^^Wv=y) 

Finally, neglecting for reasons already mentioned the squares of Tc^ and /, we find 

(31) tlo^^g^^_j^ 

(32) sm a^a^ 



I 
I 
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If wo now introduce these last values of Tio and sin a in the equation (38), we at once see 

that 

(43) Fo and consequently A^Ao 

This shows that the difference between Ti and ^o, or between F and Fo, is only of the order 
of Ic^^; that is, of the second order compared with h^, and consequently negligible. We thus 
see that in the wake of tJie aerofoil we must have nearly the same velocity as in front, hut a decreased 
-pressure. 

All the foregoing discussion brings us to the following important conclusions. Both 
assumptions for neghgible values of and/^ bring us to the same mean values of the height 
of the stream disturbed by the aerofoil and the angle a of the apparent stream deflection, 
th" variation in the magnitude of the stream velocity and stream cross section being of the 
order of A-,', and thei-efore negligible. The expressions (31) and (32), compared with the cor- 
responding expressions (22) and (23), give for the height h larger values and for the angle a 
smaller values. 

Let us now see how far the results of the foregoing discussion are verified by experiment. 

In his " NouveUes recherches sur la resistance de I'air et I'aviation," G. Eiffel, on pages 165-170, 
gives the results of measurement of the air resistance of an aerofoil disposed in the wake of 
another. Let us designate by I the aerofoil disposed directly in the wind stream, and by II the 
aerofoil disposed in the wake of I. In figure 31, Zi is the zero line of the first aerofoil and 




Fig. 31 

Zii the zero line of the second, Fo, the velocity of the stream before the first aerofoil, and i and 
ia the angles of attack of the first and second aerofoils relative to the velocity Fq. Eiffel shows 
that for the aerofoil II the air resistance is such that the angle of attack of that aerofoil, instead 
of being equal to i,, appears to be reduced to the value ir, smaller than i^. Eiffel calls i^ the 
"apparent" incidence and ir the "real" incidence. Eiffel gives the corresponding values of 
ia and ir for different angles 7 between the zero planes of the two aerofoils; and these we 
reproduce in the Table A, only referring all the angles to the corresponding zero lines. 

Table A. 



I 


■r-2* 




7-6° 




i, 


a 


i. 




a 


K 




at 


S" 

7 

9 
11 
13 
15 
17 


3' 

0 

7 

9 
11 
13 


















l',4 
2,7 
3,9 


3', 6 
4,3 
5,1 


3' 

5 

7 

9 
11 
13 


-0,4 
0,6 
1,9 
3,1 
4,5 
5,6 


3°, 4 

*,* 

5,1 
5,9 
6,5 
7,4 








3° 
5 
7 
9 
11 


-V,6 
-0,4 
0,8 
2,1 
3,2 


4°,6 
5,4 

6,2 
6,9 
7,8 

















Let us now conceive the stream velocity behind the first aerofoil as deflected and having 
the general direction F. The angle of attack, is the angle between F and Zu, so that, if a 
is the value of the stream deflection angle, we must have (see fig. 31) 

ia — = « 
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In the Table A, by aid of the values of ia and ir, the values of a are calculated for different 
values of the angles of attack i of the aerofoil I, and three different values of the angle 7, 2°, 4°, 
and 6", respectively. By using the values obtained, in figure 32 are plotted the curves oia = f{i) 
for the three values of 7 = 2°, 4°, 6°. It is easy to see that all the plotted points define well 
enough a curve, which means that behind the aerofoil I, we really have an apparent stream 
deflection depending only upon the aerofoil I, and its angle of attack, i. The angle of mutual 
inclination of the two aerofoils has no influence on that phenomenon, as ought to be expected. 
The aerofoil, II, has simply to be considered as disposed in a stream deflected by the aerofoil I, 
whose velocity in magnitude is nearly the same as in front of the aerofoil I. Knowing the rela- 
tion between a and i, we can easily calculate the air resistance of any other aerofoil, as the aerofoil 
II, disposed behind the aerofoil I. 
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i lG. 32. 



004^ 



0.03 



0.02 



00/ 



Let us now calculate according to these last data the values of the coefficient/. This cal- 
culation is made by aid of the formula (32), in the Table B, where are reproduced all the data 
concerning the aerofoil I, which are necessary for that calculation. 

Table B. 



< 


a' 


a radiants 






2 






6* 

8 
11 
14 
17 


2', 5 

4,0 

5,3 
6,3 
7,5 


0,044 
0,070 
0,093 
0, 133 
0,130 


0, 176 
0,288 
0,383 
0,475 
0,563 


0,01S2 
0,0217 
0,0308 
0,0434 
0,062 


0,00.3.88 

0,0101 

0,0178 

0,0268 

0,0366 


0,0144 

0,0116 

0.013 

0,0186 

0,0234 


0,79 

0,535 

0,42 

0,41 

0,41 



The values of the coefficient/ obtained are plotted in figure 32, where, for comparison, 
is also plotted the curve of Ic^. In tlie last column of the Table B is also calculated the ratio 
of / to fcx- These last values are very suggestive. Thej^ show that the dissipative force F 
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decreases when the angle of attack increases, and that for mean values of the last it is equal 
to nearly one-half of the drag. 

Finally I will remark that the energy dissipation, of which account was taken in the 
foregoing, and to which corresponds the dissipative force F, is that dissipation of energy 
which takes place in the immediate neighborhood of the aerofoil. It is evident that all the 
energy spent to move an aerofoil in air is dissipated in the surrounding medium; but one 
part is dissipated in the direct neighborhood of the aerofoil and corresponds to skin friction 
and turbulence connected witli it, which we evaluate by FVo', and another part is dissipated 
by the oscillatory motion of the air left behind the aerofoil— that is, by the damping through 
viscosity of the fundamental and secondary waves created. 

The experimental study of the apparent stream deflection behind the aerofoil in the 
light of the ideas here developed is very important for many 
problems connected with the design of airplanes, propeller, and 
turbines. It must, however, be remembered that we have to do 
only with an apparent stream deflection, the real motion of the 
fluid behind an aerofoil or turbine wing being generally periodical. 

Short review of some propositions on vortices. — Before 
proceeding to examine the question of the tip vortices and the 
fundamental and secondary wave, I shall state briefly some well- 
known propositions on vortices in general. 

Let us consider a small circle of radius /• lotating in its plane 
with an angular velocity w. Each part of the contour of this 
circle has a velocity equal to ru. If we now calculate for the 
contour of the circle, the quantity which in hydrodynamics is 
called the circulation — that is, the integral of the velocity v along 
the con torn- of the circle — we find (see fig. 33) 




Fig. 33. 



(44) 

or 

(45) 



2d(T 



where / is the circulation along the contour of the circle and dc the area of the small circle. 
We therefore see that the angular velocity o) of a rotating circle is equal to the circulation / 
divided by the double of the surface of the circle. 

Let us now consider a fluid element having the velocity {u, v, w) at a point (z, y, z) in a 
moving fluid mass. Let us draw through this point the axes X, Y, Z parallel to a system of 
triorthogonal immobile axes and calculate the circulation along an elementary contour with 
sides equal to dx and dy, as shown in figure 34. We have 



dl=udx + (v + ^^dx^y-(u + ^dy^dx-vdy = (^^-^^dy 



The quantity 
(46) 



^\bx byj- 



dl 
2dxdy 



which we designate by is called the component of the vortex with reference to the Z axis. 
(Compare formula (45).) In a similar manner 

, /dw dv\ , /du dw\ 

are the vortex components with reference to the X and Y axes. We have 
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The doubles of the vortex components 2wx, 2to„, 2wj are the determinants of the matrix. 



II _d d d 



II dx dy dz 
u V w 

2cOx 2co„ 2a)2 



As a consequence of continuity, the following proposition holds for vortices : 
If we have a vort«x motion at one point of a fluid, we must necessarily have vortex motion 
at all the points of a line going through that point. Such a line of small cross section do at each 
point of which the vortex has a finite value is called a vortex filament. The quantity 



is called the intensity of the vortex filament. A vortex fiilament can never begin or end in a 
fluid. It must be a closed contour or have its ends on the boundary surface of the fluid. The 
cross section of a vortex filament can be variable, but its intensity is always constant along the 
whole filament— that is to say, for a vortex filament, we have 



The vortex filament is always constituted of the same fluid particles — that is, the vortex filament 
moves with the fluid. A system of vortex filaments disposed close together form a vortex 
tube. 

Let us consider a fluid mass in motion with a vortex filament in it, and let us draw a surface 
across the fluid and take a contour on that surface. The circulation along the contour is equal 
to twice the sum of the elements of the surface multiplied by the components, along the corre- 
sponding normal to the surface, of the vortices on that surface. (See fig. 35.) 



where Wn is the vortex component normal to the considered surface. This last relation con- 
stitutes Stokes's theorem. In the application of this theorem two cases have to be distinguished. 
If the contour, by progressive shrinking, can be reduced to a point without leaving the fluid, 
the space occupied by the fluid is said to be "simply connected." If, inside the space occupied 
by the fluid, we have solid bodies or holes crossing the fluid mass, not every contour in the 
fluid can be reduced, by shrinking, to a point; and the space containing the fluid is said to be 
"not simply cormected." In a simply connected space, if the circulation along the contour 
has a finite value it means that vortices are crossing the inside of the contour and the double 
of the sxmi of the components of the vortices normal to a surface containing that contour multi- 
plied by the corresponding elements! of the surface is equal to the circulation along that con- 
tour. In a not simply connected space a finite value of the circulation can also mean that solid 
bodies or holes are crossing the inside of the contour. 

In a fluid mass in motion with vortices in it, the velocity at each point depends upon the 
distribution of the vortices. Each element of each vortex contributes to the velocity at each 
point. The components of the velocity at a point due to a vortex element ds of a vorte-x fila- 
ment of intensity / is equal to (see fig. 36^ 



where r is the distance between the point considered and the vortex element, and ip is the angle 
between ds and r. The direction of the velocity du is normal to the plane containing r and ds 
and has the sense of the rotation around ds in the sense of the vortex. The velocity at a point 
is the geometrical sum of the velocity components due to all the elements of the vortices con- 
tained in the fluid. 

"~ in the case of one single rectilinear vortex filament in an infinite fluid mass, by reason of 
symmetry the velocity is the same for all the points at the same distance from the vortex. If 



(47) 



2w d<T= I 



2co d[<7= Const. 



(48) 




(49) 
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we therefore calculate the circulation along a circle contained in a plane normal to the vortex, 
and having its center on the vortex, we get 

I = 2wd<y = 2irrv 

so that ' 
(50) v = ^^ 

That is to .~ay, the fluid velo(!ity for points around a straight vortex filament is equal to the 
intensity / of the vortex divided by times the distance r from the point considered to the 
vortex. Tlie vortex itself does not move. 

Let us consider now the case of two straight parallel and infinite vortices of equal intensity 
I rotating in inverse senses and contained in an infinite fluid mass. In figure 37 the vortices 
are represented in cross section and in plane, and L is the distance l)ctween the vortices. We 
will refer the vortices to a system of triorthogonal axes X, Y, Z as shown in figure 37, the origin 
being in the middle between the vortices and the Z axis being perpendicular to the plane con- 
taining the vortices. 

Let us first calculate the velocity at the point A, at the distance a- from the origin 0, due to 
one of the vortices. Applying the foregoing formula we have 



u 



ty-oo 



00 

dy sin <p 



and from figure (31) we easily see that we have 

\) = r cos tp; {L/2 — x) = r sin ,p 

from which follows 

dy _ —d(p 
L/2— x~ am V 

Substituting these last in the foregoing integral we get 

T 

*2 



I r 

:!2-x) / 



^ iniLI^-x) I -sin 



and integrating we finally have 



« = 4x(L/;2-x)^°"*' 



'2ir{Ll2-x) 



IT 



Considering x equal to LI2 we will find the velocity of one of the vortices produced by the other 
which is equal to 

This is the value of the velocity with which both vortices move parallel to (ho Z axis. For 
the velocity between the two vortices produced by both we will have 

(52) ,r,a — rx+i 



'2ir{L/2-xy2ir{LI2 + x) 2Tr{UI4-x') 
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For points in tlie middle between the vortices the velocity is equal to 

(53) v = ^=4u 

For other pomts between the vortices the velocity follows a parabolic law, and, as a first approx- 
imation, can be considered as nearly luiiform in the middle part between the vortices when the 
cross se-tions of the vortices arc small relative to the distance between them. (See fig. 37.) 
We therefore see that the velocity of the fluid in the middle between the vortices is exactly equal 
to four times the velocity of the vortices themselves and that the velocity of the fluid between 
the vortices relative to them is equal to twice the velocity of the vortices themselves. 
We liave now all the necessary references for the following:' 

The tip vortices.. — Let us consider an aerofoil moving with a constant and uniform 
velocity V in a fluid mass and let us designate ])y dz an element of length of the aerofoil, z 
being the distance of tlie element of the aerofoil considered from its middle cross section. (See 
fig. 38.) Let dRy be the lift corresponding to the element of the aerofoil considered . The 
quantity 

dz 

is the lift per unit of length at the distance z from the middle cross section of the aerofoil. 
According to Kutta's theorem (see Note I) we must have 

(54) ^§^ = 57/ 
so that 



(55) . /= 



hV.dz 



is the value of the circulation along such a contour as I, embracing the aerofoil at the cross 
section considered. (See fig. 38.) When z varies, and until we do not approach too near t<? 
the tips of the aerofoil, the value of / is nearly constant. As a fii-st approximation for the 
mean value of / we can take 

(56) 7=^-^ 

where L is the length of the aerofoil. Let us now move our contour to the tips of the aerofoil, 
and just before the value of I begins to change we let the contour follow the fluid in its motion. 
According to a theorem of William Tliomson (Lord Kelvm), the circulation along the contour 
moving with the fluid must be invariable, so that when the contour reaches such a position as 
II (see fig. 38) the cinulation will have the same value as in the position I. But in the position 
II we have no solid body inside the contour, and consequently we must, according to Stokes's 
theorem, have a vortex traversing that contour. 

We so come to the conclusion that vortices must necessarily run off from the tips of aerofoils, 
and that the mean value of the intensity of this vortex must be 

(57) ^=m=M^^ 

According to this statement, the tip vortices disappear (mly when the lift vanishes; that 
is to say, when the relative wind blows on the aerofoil along the zero plane. 

Having estabUshed the value of the intensity of the tip vortices, let us now consider their 
influence on the air resistance of the aerofoil. For the simplicity of the explanation we will 

> For more detailed references on vortices see the classical treatise on mechanics by Paul Appell, "Traite de mecanique ratlonneUe," Tome 
lU, Chap. XXXV, p. 389, and p. 473, § 791, and also ".\eronautics in Theory and Experiment," by W. L. Cowley and n. Levy, Chap. III. 
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put ourselves in the simplest case and will assume that the tip vortires are of small cross section 
compared with the length of the aerofoil, so that they can be considered as vortex filaments. 

If we follow an aerofoil moving with a constant velocity V, we see that the tip vortices 
run parallel to themselves, as should be the case for parallel vortices rotating in inverse senses. 
The velocity with which these two vortices displace themselves in the stream, at a distance 




Pig. 38. 



from the aerofoil where the influence of the last can already be neglected, is, according to the 
foregoing, equal to 



(58) 



and is normal to the plane containing the vortices. The motion of the fluid between the two 
vortices in the middle part relative to the vortices is equal to 



(59) 



ttL 



So that the mean value of the flow velocity between the tip vortices relative to the aerofoil is 
equal to (see fig. 39) 




Fig. 39. 

We therefore see that as a first approximation we can consider the tip vortices as bisecting the 
angle which the mean value of flow velocity between them makes with the original direction of 
the stream running on the aerofoil (see fig. 39). We also see that the angle of deflection of 
the tip vortices downward is equal to 

(60) ^""2-7 
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Substituting the value of the velocity u by its expression (58), and in the last the circulation I 
by its value (57) we find 

(61) sin^=7'^ 



If there were no other circumstances producing the deflection of the stream downward, this 
would be the correct mean value of the deflection angle. 

If now, for the calculation of the part of the drag due to the tip vortices, we should apply 
the momentum theorem, assuming that the mean value of the velocity behind the aerofoil is 
as a first approximation the resultant of the velocity V of the aerofoil and the velocity 2u = 
IlvL, we would find, according to the calculation already made on page 34, that 

a 

sm - 

.a 2_ 



and substituting in the last expression the value of sin ^ given by (61) we would get 

slcy 

where by s we have designated the quantity 
or, finally, 

which expression gives us a first approximation to the value of the drag due to the tip vortices. 
If we compare the values which the last formula gives for the drag with those values which 
direct experience shows for actual aerofoils, it will be easy to see that the drag of an aerofoil is 
much larger than that calculated by the formula (62). ' We afterwards see that the value of 
the drag given by the formula (62) decreases with the aspect ratio of the aerofoil and is equal 
to zero for an infinite value of the aspect ratio. These last conclusions are very important. 
They first show us that that part of the drag of an aerofoil which is due to the tip vortices 
practically disappears for a value of the aspect ratio greater than a certain one. We afterwards 
see that the most important part of the drag, practically the whole, for a sufiicient value of the 
aspect ratio, is due to other circumstances than the tip vortices, which circumstances are, as 
we already know, on one hand the energy dissipation in the direct neighborhood of the aerofoil — 
that is, skin friction and turbulence — and on the other, the fundamental and secondary waves. 

We can now perfectly conceive the importance of the phenomena which take place on the 
surface of contact of the solid and fluid bodies. We can not expect to be able to calculate the 
drag of an aerofoil before these phenomena are fully understood, as has been already stated in 
the preamble. 

, 

> The easiest way to see it Is the (ollowiog: 

RemarkiDg that in the formula (62) s'ky' Is a very small quantity we can write 



For actual aerofoils, as mean values we have 
So that 



6/i-l/8; — A-sO.2; r- 1/200 
2t Zi o 



Wbloh 13 oat ol proportioD with any observed value tor the drag-UIt ratio. 



\ 
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All this discussion brings us to veiy important conclusions. We have fully understood 
the influence of the aspect ratio on the drag and, so to say, the mechanism of its influence. 
When the aspect ratio increases, the part of the drag due to the tip vortices practically disap- 
pears, the influence of the tip vortices becoming negligible. The last fact taking place, the 
shape of the tip vortices is also negligible. It follows that the shape of the tips of the aerofoils 
has a negligible influence on the air resistance of the aerofoil, if only the aspect ratio has a 
sufficient value. ' 

Wo can now formulate the two following propositions: 

1. The tip vortices have an injluence on the drag of the aerofail; but this decreases with the 
increase of the aspect ratio and practically disappears for a certain value of the last. 

2. When the aspect ratio of an aerofoil has a sufficient value, the influence of the form of the 
tips of the aerofoil on the air resistance is negligible. 

These two last conclusions are in full agreement with all experience with aerofoils up to this 
time. 

We have been able, by the analysis of the tip vortices phenomenon, to undei-stand the 
reasons which require for the aerofoils a certain value of the aspect ratio and have deduced 
the slight influence of the tip forms when the aspect ratio has a sufficient value. 

Short Review of the Properties of Systems of Parallel Vortex Bows. — We will call 
vortex row a system of an infinite number of infinite rectilinear parallel and equidistant vortices 
of infinitely small cross section and of strength / equal in magnitude and sense, disposed on one 
straight line. We will give a short review of the properties of systems of parallel vortex rows, 
which play a very important r61e in the phenomenon of fluid resistance. 

The fluid mass in which the infinite parallel rectilinear vortices are considered is assumed of 
infinite dimensions in aU senses, havmg in infinity a velocity equal to zero, or moving as a whole 

with a velocity constant in magnitude and direction, inde- 
pendent of the motion which can take place inside the fluid. 

We consider the whole system of vortices cut by a plane 
normal to them and their mutual positions defined by the 
positions of their sections in that plane. 

We wLU first consider in their general outlines the condi- 
o tions which must be satisfied by a system of infinite parallel 
and rectilinear vortices in order to maintain an invariable 

configuration. 

In figure 40 are represented the sections of a system of 
vortices, as above mentioned, cut by a plane normal to them. 
The vortices are assimied to maiiitain for one moment an invariable configuration and are 
referred to the orthogonal {X, Y) axes moving with the vortex system. Let us concentrate 
our attention on one of the vortices of the system, say A, with the coordinates x and y, and 
consider this vortex A undergoing an infinitely small displacement Sx and By relative to the 
other vortices, the last maintaining their configuration, and let u and be the components of the 
velocity of the vortex A. As is well known, u and v are the components of the velocity pro- 
duced at the point x, y by the other vortices of the system. The displacement of the vortex 
A will produce the variations du and dv of its velocity components, which will represent the 
components of the velocity of the vortex A relative to the other vortices, and we will have : 

, bu . bu . 
ox oy ^ 

dv ^ bv , 
bx by •' 

The vortex A wUl be stable for every virtual displacement (5a;, dy) provided du and dv, respec- 

i And i( the tips do not have any extravagant form which could give rise to additional local aix resistance. 
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tively have signs opposite to the signs of the displacements bx and by. But it is easy to see 
that, so long as 

bu bu bv bv 
bx' by' bx' by 

have finite values, it is always possible to imagine a displacement (5a;, by) for which du and 
dv will have the same signs as bx and 5?/. In other words, a stable displacement of a single 
rectilinear vortex among others maintaining their configuration appears as impossible. 
But in the case of 

(63) ^ = 0;t^ = 0:|^ = 0;^ = 0 

^ ' bx ' by bx ' by 

we have 

du = 0; dv = 0 

for all infinitely small virtual displacements of a single vortex. Under such conditions the 
vortex considered will be in a state of neutral relative equilibrium among the other vortices. 

As for the space between the vortices, we must have the equation of continuity satisfied; 
that is, 

bu bu Q 
bx by ' 

as well as the vortex intensity equal to zero 

bu bv _Q 
by bx~ ' 

The four conditions (63) are reduced to two conditions. It is thus sufficient to have two of 
the quantities (63) equal to zero in order to have a neutral state of equilibrium of a single vortex. 

It is easy to see that the conditions (63) mean that the velocity (u, v) must have a maximum 
or minimum at the point where the vortex A is disposed. 

When the relations (1) are satisfied, a small displacement of the vortex A does not produce 
a change of its velocity. But let us consider the influence that the displacement of the vortex A 
can have on the other vortices of the system, say on the vortex B, for example. Let us imagine, 
first, that the vortex A undergoes a displacement Sr along AB only (see Fig. 41 ). Let 7 be the 
strength of the vortex A. Before its displacement the vortex A was producing in 5 a velocity 
normal to AB, that is, to r and equal to I/2irr, which, combined with- the velocities that the 
other vortices of the system produce in B, keep this vortex in relative equilibrium. After its 
displacement, the displaced vortex A will produce in B a, velocity equal to //2ir (r + br). We 
thus see that we can consider the displacement of the vortex A as producing in B an additional 
velocity equal to 

^^^^ 2Tr(r + br)~2^^~2^^'' 

which is normal to r (see Fig. 41). In the same manner it will be easily seen that the displace- 
ment 571 of A normal to r will produce in B an additional velocity directed along r and equal to 
(see Fig. 42) 

(65) -~^bn 

If we now let B displace itself in the direction of the additional velocity communicated by 
the displacement of -4, we will at once see that such a displacement of B will produce in A the 

additional velocities — along r — ^Sn along the normal to r where ic is a factor of pro- 

portionality, which will be exactly inverse to the displacements br and bn when the two vortices 
considered are of inverse sense (see Fig. 43), and which will have the same senses as 5r and bn 
when the two vortices considered have the same senses (see Fig. 44). 
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The foregoing preliminary discussion is given here only to show the nature of the question. 
It will, however, allow us to arrive at a general understanding of tlie stability conditions of 
vortex rows. 

Let us consider first one vortex row, i. e., all the vortices in the same sense (see Fig. 45). 
If we disj)lace one of the vortices of the row, say A, fii-st to the position A', considering the 
other vortices for one moment immobile, it will be easy to see that in the position A' this vortex 
will get from the other vortices an additional velocity in the sense oi the arrow 1 . If we displace 
this same vortex to the position A" it will get an additional velocity in the sense of the arrow 2. 
Thus for neither displacement will there appear an additional velocity directed toward the 
original position of A. If we now consider the additional velocities whicli the vortex will I'eceive 
from the displacements of the other vortices produced by the displacement of the vortex A, 
it will be easily seen from the foregoing, since all the vortices have the same senses, that these 
additional velocities will only increase the original displacement of the vortex A. Thus one 
vortex row appears as an unstable configuration. 

Let us now consider two parallel vortex rows rotating oppositely in each row. The 
additional velocity given each displaced vortex will consist of the velocity due to its own dis- 
placement and of the velocities this vortex gets due to the displacements of the other vortices 
produced by its own displacement. This last additional velocity which is caused by the 
vortices of the same row as the vortex considered, has a destabilizing action, as all the vortices 
of the same row rotate in one sense, but the additional velocity caused by the vortices of the 
other row will produce a stabilizing action, as they rotate in inverse sense. We thus can conceive 



d\ / -A ^^'^r^ 

Fig. 47.' 

that two parallel rows of vortices with inverse rotation in each row can have a stable configura- 
tion. The investigations of Karman have shown that two parallel rows such as represented in 
Fig. 46 can not be stable, but two rows in quincunx as represented in Fig. 47 can be stable for 
a certain value of the ratio d/2l—\ of the distance d between the two rows to the distance 21 
between the vortices in each row. The most probable value of X for a stable configuration of 
two vortex rows in quincunx seems to be ' 

(66) ^^ = X = 0.283 

If a single vortex row could be stable it is easy to see that it would be immobile, because 
the velocities which all the vortices communicate to one of them mutually balance. As for 
a quincunx arrangement of vorte.x rows, it is easy to see that it will move with a constant 
velocit}' parallel to the general direction of the rows (sec Fig. 47) because each vortex of one 
row will receive no velocity from the vortices of the same row, but from each two vortices of 
the other row, disposed symmetrically relative to the vorte.x considered, it will receive a resultant 
velocity directed along the direction of the row. 

Let us calculate the value of this velocity. For that purpose let us first calculate the 
velocity produced by a single vortex row in a point P. Let us refer the vortex row 

A's, A' I, Ai, A3, A^ to a system of orthogonal axes X'O' Y' whose X' axis is parallel 

to the vortex row considered, and let us also consider the system of a.xes XO Y parallel to the 

I V. Karman has published two papers on that question. In his flrst paper "Uber den Mechanismus des Widerstandes, den ein bewegter 
Korper in elner Flussigkcit erfahrt" Nachrichten von der Konlglichen Gcsellschaft der Wissenschaften z\i Oottingen, 19U, p. 509, Karman found for 
X the value 0.367. N. Joukowski has made some experiments which seem to verify this value. See his " Aerodynamique," Paris, 1918. p. 20'.. In 
his second paper "Uber den Mechanismus des Flussigkeits und I/uftwiderstandes" Physikaliseho Zeitschrift, January 15, 1912, Karman flnds for 
X the value 0.283 and produces the results of experiments, also verifying this last value. This question thus demands further investigations. 
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system X'O' Y' but whose origin 0 is midway ])etwecn the vortices A\ and A (sec fig. 48). Let 
ij) ho the coordinates of the ])oint P in reference to the axes X'O' Y' and (x, y) the coor- 
dinates of the same point in reference to the XO Y axes. The coordinates of the origin 0 in 
reference to the axes X'O' Y' will be designated by r^„, jjj). We have 



(67) 



The coordinates of any vortex A of our row in reference to the X'O' Y' axes will be designated 
by (a, b). As is well known, we have ' 



(68) 
with 



u- 



.^^ I 1 



z = ^ + iri; a=a + ib;i= V"^ 



and where u and v are the components of the velocity produced in the point P by the vortices 
of the row considered, / the intensity of each vortex of the row. 



._o- 



— o, 



-<7 



\V 

!/ 

i/ 

1/ 



.io -p- 



Fig. 48. 



For the vortices of our row disposed symmetrically relative to 0, the quantity a has for values 



a3=(lo + 3Z)+i'7o 
«6=(^o + 50+ii?o 



a 
a' 

a', 



(f„-3Z)+ii7o 
(?o-5Z)+i7;o 



Thus 

Z-a, = (2 — i)+t1/ 2-a',= (x + O+i^ 

2 — ^3= (x — 3Z) 4-12/ 2-a'3= (x + SZ) 

z-ai={x-5l)+iy z-a'^={x + ol)+iy 



and, consequently, 
with 

we thus see that 



1 _ 2(.x+ iy) _ 2p 



2— a, 2 — a'i {x + iyY — P 



p-x + iy 



2p 



_ p='-(2A:-l)2P 



1 See for example "Traitc de Mecanique Ratlonelle," by Paul AppeU, Vol. HI, p. 481, 
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If we now remember the well-known formula 

k = ao 

. . yi 4-20 

k=l 

and identify this formula with 

k-oo k-K A 



k=i — (2i:-l)V-^ 

it will be easily seen that with 

^-21 



we get 

k=oo 

2Z 



and thus 

I . vp 

or 

IT . . IT 



l-2tg^ a;tg|^% 



-m-« = 4jjtg 22=5^tg-^^ 



^ / tg^^x+itgh^^y tg x-+Ughy' 

'S-2Ugf^xtgh|^, 4Zl-2Ug.'tgh,' 



introducing the notations 

"'^21*' ^ ^2Z 



, IT / T 



Thus 



_ . _ _ 7^ tg z'+t tgfa _ tga;'(l-2 tgh'y') . tghy'(l +2tgV) 
^ ~ 4Z 1 - 2i tg a;' tgh y ' ~ 1+ 4 tg'x' tgh'y ' * 1 + 4 tg^x' tgh^y ' 



and we finally find 

rfio, / tgx'(2tghY-i) / 21 [^*g^' 5-0 

^ ^ ^ 1+4 tg2 22 (^-fo) tgh^ 2i^''-''»^ 

7 tghy-(l+2tg'x-) / fcg'^lz(''-''o)[^+2^g'|] (^-^o)] 
4Z l+4tg='x'tgh^y' 4Z , , . ^ , x ,^ > ^ . uj , . 

These last formulae constitute the general expressions for the components u and v of the velocity 
at a point P tj) produced by a single vortex row. 
I will here note that, since for 

(71) t9hHr,-Vo)^i= 1/2. 

we have 

v = 0 

for any value of |, a single vortex row has around it such a flow that at the distance (it — Vo) 
fixed by the relation (71), we have two streamlines that are straight lines. 



\ 
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If we now apply the formulae (69) and (70) to a system of two vortex rows in quincimx in 
order to find the velocity communicated to one row by the other, as we have i = {o for each 
vortex, we find 

u = 0 

and ^ 

(72) 2?^''"''"^^ ~4Z 21*^ 
where 

v-Vo = d 

and for the stable configuration of two vortex rows in quincunx with X = (Z/2i = 0.283 we find 
for the magnitude of u the value 

(73) usO.35 ^ 




Fig. 49. 



Let us further calculate the momentum q counted per unit of vortex length that corresponds 
to two vortices A and B of equal intensitj^ / rotating in inverse senses which we wUl refer to 
the system of XO Y axes. (See Fig. 49.) It is easy to see that the component of the momentum 
of the two vortices in the direction AB is equal to zero. Because if we consider a fluid strip 
parallel to the direction of AB such as the component of the momentum along the du-ection 
of AB of the fluid element situated at a point such as is equal and directly opposite to the 
momentum of the fluid clement p'l symmetrical to pi in reference to the axis Y (see Fig. 49). 
The component of the momcntiun along the direction of AB which corresponds to the whole 
strip hfii is thus equal to zero. As the same takes place for any fluid strip parallel to AB, 
the resultant momentum corresponding to the two vortices A and B will be normal to AB. 
For the calculation of this momentum, let us divide the fluid into strips normal to AB such as 
6j6,. The component of the momentum of the fluid element situated at a point such as 
along the normal to AB, and due to one of the vortices is equal to (see Fig. 49) 

. , , / m-x 



I 

r 
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and for the whole strip considered 



when we note that 
We thus find 




FlQ. 60. 

which quantity appears to be independent of the coordinate x. It will be easy to see that this 
hist momentum is positive for the vortex B for all the strips at the left of B and negative for all 
the strips at the right of B. The inverse takes place for the vortex A (sec Fig. 50). Thus when 
we calculate the total momentiun g, the momentimi of the strips outside A and B will mutually 
cancel, and there will be left the momentum of the strips between A and B. We thus find 

(74) 2 = 2 J • h = 8lh 

where h is the distance between the two vortices considered. 

Let us now consider a system of two parallel rows of vortices in quincunx. We can always 
conceive this system as built up by the superposition of two identical vortex row systems with 



I 

i 

I 



54 



ANNUAL REPORT NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS. 



the intensity of each vortex equal to 7/2. (See Fig. 51.) To each pair of vortices of intensity 
7/2 will correspond a momentum normal to the line joining them and equal to 

q = bli ^ 

If wo consider the sum of all this momentum for the whole system, it is easy to see that we will 
get only a resultant along the general direction of the rows. Each pair of vortices of intensity 
7/2 will contribute to this resultant by a momentum equal to 

(75) . «aM=«4 

and the resultant momentum counted per vortex of intensity 7 and per unit of length of the last 
wUl thus be equal to 

^^2 




Fig. 51. 

Summing up from the foregoing all the data relating to a system of two parallel vortex 
rows in quincunx we see that — 

The quincunx system maintains a stable configuration for 

(76) 1^0.283 

The system communicates to itself a velocity parallel to the general direction of the rows 
equal in magnitude to 

(77) u = ^^tghg^0.35^ 

The resultant momentum of the system is directed along the general direction of the rows 
and counted per vortex and per unit of length of the last, is equal to 

(78) hd I 

The Fundamental and Secondary Waves. — We will now make an attempt to calcu- 
late the order of magnitude of the fundamental and secondary waves. As we have seen in the 
foregoing, the fundamental wave will be produced by the vortices in quincunx built from the sur- 
faces of discontinuity which are, in the limiting case, the boundaries of the wake behind a body 
or aerofoil. The secondary wave is produced hy the vortices in quincunx built from the 
KirchhofT-Lord Rayleigh surfaces of discontinuity. In some cases, in all probability only one 
of the waves will be formed; in others both waves will appear simultaneously and one will 
propagate itself in the other. The study of the conditions of formation of the fundamental 
and secondary waves demands further deep investigations. I will here consider only the cases 
in which each kind of wave appears separately. 
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Let us consider a, solid body or an aerofoil disposed in a uniform fluid stream having a 
general velocity equal to V. If a system of two parallel vortex rows in quincunx is assumed 
to appear behind, the momentum counted per unit of time, corresponding to the vortices 
appearing, will be equal to the momentum per vortex, multiplied by the number of vortices 
formed. As the vortex in the flow will have a velocity equal to u the number of vortices 
formed will bo equal to 

and the corresponding momentum equal to 

(79) ^^^W^- 

To a first approximation we can assume the drag of the body or aerofoil considered (counted 
per unit of length) equal to the last quantity — that is ' 

(80) R^ = SdI^~^^Jc^SbV' 

or, introducing for u and d their values in (77) and (76), we get 

(81) Av6 = 0.283 -^(^1-^^) 

The following considerations allow us to estimate the value that the intensity of the vortex 
in quincunx built behind a body or aerofoil may have. 




Fig. 52. 

Let us consider in a fluid flow a surface of discontinuity, or vortex sheet, on both sides of 
which we have a finite velocity difference equal to w' = V^ — (sec fig. 52). For a contour 
such as C drawn between two points A and B of the surface, whose distance is 21, the circulation 
will have the value 

2lV,-2lV^ = 2l{V,- V^)=2lw' 

When now the vortex sheet considered goes over into a row of vortices whose mutual distance 
is 21, in the ideal case the circulation will remain invariable, and if a vortex is built between the 
points A and B, the intensity of it will be equal to 

(82) I=2lw' 
Thus in the case of the fundamental wave we will have 

(83) w'=w and I=2lw 

where w is the mean wake velocity; and in the case of the secondary wave we will have 

(84) w'= Fand /=2tF 

the fluid inside the Kirchoff-Lord Rayleigh surfaces of discontinuity being at rest and the 
velocity on the surface being equal to F. 

» Tlifi whole purpose of this paper is to be first of all quite elemciU.iry. That is why I have allowed myself to give tlie forej:oinR formula, 
which is not quite exact, tlie flow periodicity behind the body and the pressure dillcrencc in front and behind ha\dng been neglected. The error 
committed usins this formula can reach 30 per cent; that is why this formula Kives only an idea of the order of magnitude of the quantities considered, 
which is the only thing we wish to reach here. 
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Indoducing the value (82) of / in (81) we find 

M = 0.283 • 21 y (^1-0,35'^) 
and thus for the wave length 21 we get 
(85) 21- 



0,283y (1-0,35 



(86) 



Thus at first approximation the length of the fundamental wave appears equal to 



21- 



0,283 p(l-0.35p) 



and the length of the secondary wave equal to 

(87) 2Z^5,4 k^b 

To give a concrete example, let us take an aerofoil with 

6 = 2 mi.; tj. = 0.05 

and assume the ratio w/V to be of the order of 0,01.' 

In such a case, for the length of the fundamental wave we find 

07 0.01 „^ , 



and thus 



and 



0.283 • 0.01 

(Z = 0.283 • 2Z^10 mi. 
For the length of the secondary wave we get 

2Z = 0.1 • 5.4^0.5 mi 

(Z^O.15 mt. 



Secondary waves have been observed by several experimentators, and are fully of the 
order of size given by the foregoing formulae. But fundamental waves, as far as the author 
has knowledge of the subject, have never been experimentally observed. The scale of this 
phenomenon shows the great interest of its experimental study. The possibility of the for- 
mation of the fundamental waves explains the action which bodies moving in a fluid may 
have on each other when (hey approach one another, as has been observed in some cases be- 
tween airplanes and sliips. The phenomenon of the fundamental wave indicates also how 
complicated is really the comparison of the fluid resistance of a body moving in a free and 
in a closed .space. I will, finally, mention once more (he fact that the phenomenon of the 
fimdamental wave is a consequence of the fluid viscosity and can not be conceived in an ideal 
fluid. 

The Pressure Distribution on the Surface of the AerofoU. — One question has 
been left so far without discussion: It is the pressure distribution on the surface of the aerofoil. 
The general outlines of this phenomenon are easy to underetand. The velocity of the flow 
running on the aei'ofoil is increased above the aerofoil and decreased below, which has as a 
direct consecjuence the decrease of pressure on the upper surface and the increase of pressure 
on the lower surface of the aerofoil. But we are not able to make the exact calculation of the 
pressure distribution along the surface of (he aerofoil. The pressure distribution is very 
closely connected with the phenomena which take place on the surface of contact of the fluid 
and solid body. Until some new conception throws a new light on these phenomena, it does 
not seem that the pressure calcula(ion can be started with am^ success. The general ideas 
which were developecf in the preamble indicate the way in which the solution of the problem 
will probably be found. 



' The value ottcIV is connected with the value ot kz. The elementary considerations developed in the beginning o( this chapter may be used 
to and the mutual order of magnitude of kx and vol V. 



NOTE I.— Generalizalion and General Discussion of Kutta's Theorem on Circulation. 



The circulation thoorem which 1 have in view in the present note was (irst indicated for a particular case by W. M. 
Kutta." Soon afterwards, Kutta and Joukowski ^ have recognized the generality of this theorem. This theorem 
is announced as follows: 

When a rectilinear and uniform fluid current, having at infinity the velocity T, flows normally to the generatrix of an 
infinite cylinder from any section, and when the circulation along a contour embracing the cylinder and situated in the plane 
of one of its orthogonal sections has a finite value I. the component Hy of the resultant force of the fluid on the cylinder, taken 
along the normal to the velocity and referred to a unit length of the last, is equal to the product of the velocity V, the circulation 
land the density of the fluid: The sense from Ryto Vis coincident with the sense of the circulation. 

According to this theorem, the lift produced by a unit length of the cylinder considered is expressed in magnitude 
by the following formula: 

Ry=SVI 

We shall establish two fundamental and quite general relations from which the circulation theorem will appear as a 
particular case. 



y 




X 



FlQ. X. 



di/>=udy —vdx=:^^y-^^dx 



Let us embrace the infinite cylinder considered by any contour disposed in the plane of one of its orthogonal 
sections. Let Wbe the velocity of the fluid at the point (x y) of the contour; u and v the components of the velocity 
TF along the axes (see fig. 1); dx and dy the projections of the element of the contour on the axes. Let us designate by 
X and Y the components of the resultant force of all the exterior forces applied to the fluid contained in the contour 
considered and let us apply the theorem of momentum to the motion of the portion of the fluid considered. We then 
have 

(1) Y=fvdm; X=fudm, 

the integral being taken around the contour and dm representing the fluid mass which flows out per unit of time through 
an element of the contour. Let us designate by <l> the current function. By the definition of that function, we have 

(2) dm,=idyp 
and also 
(3) 

w 

with 

i>y ' dx 

Substituting in the first of the equations (1) the value of dm taken from the equation (2) we obtain 

(5) Y=fSvd\f>=fSv{udy-vdx) 
or, 

Y=fSli:(udy-vdx)+u^dx-u^] 
=fiu{udx+vdy) —jS(u''-{-v'')dx; 

and, remarking that 

(6) udx+vdy=dl 
I \V. M. Kutta. lUustrirle Aeronautische MitleilunKen, 1902. 

' W. M. Kutta. Sitzungsberichte dor Konislidien Baycrisdien .Vkademie der Wissenschaflcn, Miiiiclicn, 1910 and 1911. 
' N. E. Joukowski. tieometriche Untersuchungen Uber die Kutta'sclie stromung, Moscow, 1910, 1911. See also his coarse " Aerodynamique,'' 
Paris, 1916, p. 139. 
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is nothing else than the flow dl along an element of our contour, in a counterclockwise direction, i. e., a direction such 
as will turn the axis of X into that of F, we get 

(7) Y=^iudl-^m'dx 

and tiually, iut^rating by parts the first term of the second member of that relation, we get 

(8) Y= [«u /- J«/ (/«] - Js W^x 

which expression holds for any contour and constitutes the first of the relations we wislied t>( get . 
Applying that relation (8) to a contour along which 

v=o; u= V=€onst. 

we easily see that we have 

(9) jildu=0: ^bWMx=0 
and consequently T reduces to 

(10) Y=i VI 
I being the circulation along the contour considered. 

Following the same way with the second of the equations (1), we get 

(11) X=jsud<;,=jSu{udy-vdx) 

= J J (u'dy —uvdx—v'' dy+v'dy) 
=j5(u^+v'')dy- jsv{udx+vdy) 
<12) X=fsWdy-jsvdI 

(13) X=js WHy- /- J6 /</«] 

The last of these equations constitues the second relation we wished. 
Applying this last relation to a contour along which 

,.=0; u= l^=Const. 

we easily see that we have 

(14) X=0; 

all the three termsof the second member of the relation (13) being e()Ual to zero. 

Let us now stop to note the exact interpretation of the relations (10) and (14). .\s it has been indicated. A' and }' 
are the components of the resultant forces of all the exterior forces applied to the volume of fluid contained in the con- 
tour considered. These forces are: First, the pressures of the cylinder on the fluid, which are equal and opposite to the 
pressures of the fluid on the cylinder; second, the exterior pressures on the contour. Let us consider a contour over 
which v=o; u= K=Const. and which is limited in one sense by two stream lines sufficiently distant from the cylinder 
so that they are parallel to the X axis, and in the other sense by two lines perpendicular to these stream lines. Along 
the stream lines parallel to the X axis we can consider the Bernouilli constant as being effectively constant, and in con- 
sequence the pressure p constant and eijual to the exterior pressure p„, the velocity Fbeing constant. Under this condi- 
tion the component along the Faxis from the exterior pressures on our contour will be zero, and F will represent the 
negative of the component of the pres!3ures of the fluid on the immersed cylinder. The expression (10) consequently 
gives us for the numerical value Ry of the lift of the fluid on our cylinder R^—b VI. But if we consider a stream line 
which flows near our cylinder, there must be some interior lo.ss through viscosity along this stream line because each 
immersed body gives rise to drag. The Bernouilli constant along such a stream line must necessarily decrease, and 
when we reach the side of the contour, parallel to the l^axis, where the velocity 1' has again taken its original value, 
the pressure there will not take its original value the Bernouilli constant having decreased. The relation (14) 
consequently expresses the fact that the component along the X axis of the resultant force of the exterior pressures on 
our contour is exactly equal to the drag, and this still in the case when the sides of Ojur contour are moved to infinity. 
In the last case, the exterior pressures tend to their limit value p„, but this is not reached ; and the integral 

always remains exactly equal to the drag. Messrs, Kutta and Joukowski, who were the first to establish the relations 
(10) and (14), have limited themselves to a consideration of a perfect fluid. In that case, having no interior losses, the 
Bernouilli constant has an invariable value along any stream line and the relation (14) expresses, then, the fact that 
the drag of an immersed cylinder is zero. But it is absolutely unnecessary to limit ourselves to a perfect fluid, since 
the theorem of momentum from which the equations (10) and (14) are a direct consequence, is applicable whatever 
the interior forces acting on the considered system are. 

We are thus brought to the general conclusion that, for any contour surrounding an immersed cylinder, the fol- 
lowing general relations must hold : 

(15) jpdx-R,= jsv{udy-vdx)=js udI-JsWVx=[Su I-js Idu]-jsW'dx 

(16) ^pdy-R^=jsu{Tidy- vdt) = Js W'dy -jsid I=Js Wdy - [5 !- I-js Idv] 

which connect the lift and drag of the cylinder, referred to a unit length of the \ast, with the flow around this cylinder. 
In the application of these formulas, three particular cases have to be distinguished: 



I 

i 
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iBt. The formulas are applied to the contour of the cylinder itself. The contour of the cylinder being a stream 
line through which we have no flow, we must simply have 



which is the case considered in classical hydrodynamics. 

2d. The formulas arc applied to a contour wliich consists of stream linos and normal lines. (For the definition 
of the last lines, see Note II.) In that case the integrals which figure in the second membersof the relations (to) and 
(IG) have to be calculated only along the normal lines. 

3d. The Kutta case 





NOTE II. — Generalization of the Bernouilli Theorem. 



For the determination of the pressures in a fluid, wo have tlie liernouilli theorem, which furnishes us the law of 
variation of pressure along a stream line and also along a vortex line. Wc also know that the Bernouilli theorem is 
applicable to the whole fluid, considering the Bernouilli constant as invariable when tho fluid motion is irrotational . 
But in the general case, when we go from one stream line to another, the Bernouilli constant changes its value. What 
is the law of variation of the Bernouilli constant in the whole fluid mass in the general case? It is the solution of this 
important question which the present note gives. We so arrive at the general solution of the problem of the distribu- 
tion of pressures in a fluid mass. 

Let us consider a fluid mass in a steady state of motion. Let us consider in this fluid mass the stream line curves 
and also two other families of fundamental curves: the normal lines, defined by the property that the tangent at each 
point to those curves coincides with the piincipal normal of the stream line passing through this point, and the binomial 
lines, defined by the property that the tangent at each point coincides with the binormal to the corresponding stream 
line. The stream lines, the normal lines, and the binormal lines form a system of triorthogonal curves. 

These curves have for equations: 

The stream lines: 

u V w 

The normal lines: 

dx dy _ dz 

Bxv—Cv~Cu--Aw Av — Su 

The binormal lines: 

dx_dy_dz 
A~ B~C 

In these equations u, v.warethecomponentsof the velocity of the fluid and C the determinants of the matrix: 



u 


V 


w 


du 


dv 


dw 


dt 


dt 


dt 



For example 

dw dv 

expressions in which j". ^ are the total derivatives; for example, 

du du dv dw 

the motion being steady. 

Let us consider a fluid element contained in the elementary parallelepiped, whose edges dr, du, df} are respectively 
directed along the stream lines, the normal lines and the binormal lines. On these curves, we get the following positive 
senses: On the stream lines, the sense of the velocity of the fluid particles; on the normal lines, the sense toward the 
center of curvature of the corresponding stream lines; on the binormal lines, the positive sense is chosen in such a way 
that the trircctangular trihedral (dr, dv, d0) be positive. 

Let us apply the d '.\lembert principle to the fluid element dr, dv, dff and let us consider for the sake of simplicity 
the fluid as incompressible and having no weight (see fig. 1) . The resultant of the exterior pressure on the fluid element 
has for components: 

— ^ dr dv dp along dr 
Or 

-^drdvdff -. - dv 

Ov 

-^J'drdvdp -.-dff 
Ov 

p being the pressure at the point considered. 

Tho resultant of the forces of inertia applied to that element has for components: 

— '^Sdr dvdff along dr 
at 

— Y^Sdr dvd0 -. - dv 
p 

0 -. -dff 

i being the density of the fluid at the point considered; Vthe velocity, and p the radius of the principal curvature. 
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(1) 

(2) 
(3) 



According to the d'Alembert principle, we must have: 



= 0 



Oi/ p 

dp 

This system of relations represents the eq uation of motion of the fl uid referred to the triorthogonal curvilinear system o f 
stream lines, normal lines and binormal lines, which can l)c called the natural cunilinear coordinates of the fluid or. 
shorter, the natural coordinates of the fluid. 

I. The equation (1) brings us directly to the Bemouilli theorem. We have 

dp dV dp dVdr ^P,.^dV 
d^+^'di^d'r + ^dFdt=d~r+^^d^ = 0 

and integrating along a stream line, we get 

a relation which constitutes the Bemouilli theorem, C being the Bemouilli constant. 




/ 
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/ / 

// ' 
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— 1 
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1 X 



Fio. 1. 



FlO. 2. 



II. The equation (2) gives us the distribution of pressure along the normal lines. Int^;ratiag this equation along 
the normal lines, we get 

(5) 



Jy2 
i—dv=C 



This last equation is susceptible of the following important transformation: 

Let us designate by u,, co^ the components of the vortex and by F^, F,, Fp the components of the resultant 
velocity F along the directions dr, dv, dp at the point considered. We have 

Vr=V; V,=0; Vb=0 

The relations between the double of the component of the vortex u> and the component of the velocity are given by 
the determinants of the matrix 





d 


d 


d 


(6) 


dr 


dy 













We thus have 

(7) 
or 

(8) 

On the other hand (see fig. 2). 
(9) 



„ ^F, dVr 



dF, 

dr 



-2w, 



dr 



jyvrjyv 

dv dv 

Vde^v 

dr p 

,dVy . 



do being the contingence angle. Substituting this last value of -g^ in the relation (8), for an integration along a normal 
line, we get 

(10) dv= - 



dV 
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and substituting this last value of di/ in the equation (5), we get: 

J SVdV 

Finally, the integral of this last relation is susceptible of the following transformation 

r sVdv _ pP-^^[i-(i-2pf)+(i-2p^0] 

Jl-2p^^ J l-2p"/ 

/r,VdV-.2^ 
svdV- I 

l_ r svdV 

J dVdV 
2pti)8 



2 



and the equation (11) takes the form 

(12) p+-2-=C+ 

2pta0 

which relation gives the distribution of pressure along the normal lines. 

We easily see that the last equat ion has the form of the Bernouilli equation, only the integral which figures in the second 
member determines the variation of the Bernouilli constant when we go from one stream line to another along a normal line. 

If we put 

(13) Ac= 

J 



2pu0 ^ 

the equation (12) takes the form 

(14) p+*I^=C+AC 

We now see, it is suflicient for w/s to equal zero along a normal line — which means that on the normal line consid- 
ered the vortex w is in the contingence plane — for the integral AC to be equal to zero and therefore for the Bernouilli 
constant to be invariable along the normal line considered. It is evident that is zero when a)=0. 

The integral AC can be written in the form 

J idV 
i 
2po)/3 V 

and is then susceptible of the following geometrical interpretation: The denominator of this integral represents the 
difference between the inverse of the speed which the fluid particle would have if rotating with the angular velocity 
2u/9 around the center of curvature of its instantaneous position and the inverse of the resulting velocity V of the 
particle. 

III. The integration of the equation (3) along the binormal lines brings us directly to the conclusion that along 
those lines 

(16) p=Const. 

that is to say, in the case of a nonheavy fluid, the binormal lines are isobars. It will be easily seen that in the case of a 
heavy fluid, tlie distribution of pressure along the binormal lines will be the same as if the fluid were immobile. 

We also see that for the case of irrolational motion of a fluid the binormal lines are also the lines of constant velocity, 
the Bernouilli theorem being applicable to the whole fluid mass. 

The system of relations for (11), (12), and (16) fully determines in the general case the distribution of pressures in 
a fluid mass in motion. This system of relations carries us to the following important consequences, which I will 
indicate in general outlines: 

I. It is sufficient to know the distribution of pressure along a surface cutting all the binormal lines in order to know the 
distribution of pressure in the whole fluid mass. 

This proposition is a direct consequence of the fact that the pressure is constant along a binormal line. 

//. On both sides of a vortex layer, even thin, there can exist a difference of pressure which can be of sensible value. 

To convince ourselves of such a possibility, it is enough to picture a vortex layer in which the quantity 



V-2poiB 



AN INTRODUCTION TO THE LiVWS OF AIR RESISTANCE OF AEROFOILS. 



63 



has a small value inside the layer, which can happen without V and we having excessive values. Then, when trav- 
ersing the layer along a normal line, the integral 

- r SVdN 
J V-2pwe 

will have a large value, and, consequently, according to formula (11), the difference of pressure on the two sides of 
the layer can have a sensible value. 

The conception of a thin vortex layer maintaining a sensible difference of pressure is very important for the under- 
standing of many hydrodyoamical phenomena which take place in real fluids. I shall take a typical example. 




Fig. 3. 



Let us consider a propeller having any number of blades and working in free air at a fixed point, for example. 
"uhe ntation of the propeller creates a well-enough limited fluid stream. Let us follow a stream line in the sense of 
/he m)tion of the fluid projected by the propeller (see fig. 3). When we reach a point such as B disposed before the 
propel|fr, the pressure p must be necessarily less than the exterior pressure Poi because the velocity is all the time 
increasing when we approach the propeller and at points such as J., and we have pressures very close to the 
pressure p. But when we go through the plan of rotation of the screw, the pressiu-c increases and in a point such as 
C disposed directly behind the propeller, the pressure p' is generally greater than the exterior pressure p^. It 
would 1)0 difficult to conceive the existence of different pressures p' and p at points C and A.,, if it were not for the 
vortex layer, which consequently must constitute the surface of the fluid stream created by the screw and which is 
capable of maintaining differences of pressure. Without the knowledge of the existence of the vortex layer con- 
stituting the surface of the stream created by the propeller the distribution of pressure around the propeller would 
be difBcult to conceive. 



NOTE HI.— The Equation of the Metacentric Curve. 



The aerofoil considered is referred to any system of AT and Faxesinvariablyconnected with the aerofoil, for ex- 
ample, the one represented on Fig. 1. The air resistance R of the aerofoil is resolved into two components R'^x and 
fi'j, along the X and Faxes. These components are connected with the drag Rx and lift Ry of the aerofoil by the 
realations 

(1) R'x=Rx cos a—Ry sin a 

(2) R^y^Rx sin a+Ry cos a 

Let us designate by N the moment of R relative to the origin 0. We have 

N=R. I. 



(3) 



The quantities 



R'x; li'y: ^ 



fully define the vector R in position and direction and are sometimes called the Plucker's coordinates of a vector. In 
these coordinates the equation of the direction of R is 



(4) 



X R'y-y R'x-N=0 




\ 



In fact, the normal form of the equation of the direction of R is 

X coe ip+y sin <p—C=0 i 



but 

and thus 
consequently 



R^y=R coe <f>; R'x=—R sin >p 

R'v ■ R\ 

coBip=--j^; sm <p= — ^ 

* ^ R 



X R'y-y R^x-Rl=0 

and the equation (4) is thus established. 

The metacentric curve is the envelope of the consecutive positions of the air resistance R^ and thus is fixed by 
the system of relations 

(5) I xR^y-y R\-N=0 

da 



(6) 



■ y dR'x ^ 



da dot 
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from which we find 



TV 


-R\ 




R\ 


N 


ON 


dR'x 




dR% 


dN 


da 


da 


(9) y= 


da 


~d^ 


«v 


-R'x 


R\ 


- R\ 


dR'p 


dR'x 




dR'y 


dR'x 


da 


da 




da 


da 



which represent the equation of the metacentric curve iu parametric form, a being the parameter 
These equations can also be written 

^^~da daWxJ 

"sa da jAr'vL - _mJjvi'lrvJ 

'^^'da-^y-d^ OaKR'yJ daKR'yJ 

which can be used for the tracing of metacentric curves. 
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NOTE IV.'— The Mechanism of Fluid Resistance.' 



By Th. v. Karman and II. Rubach. 

The resistance of a solid body moving with a uniform velocity in an unlimited fluid can be calculated theoretically 
only in the limiting cases of very slow motion of small bodies or of very high fluid viscosity. We are brought in such 
cases to a resistance proportional to the first power of the velocity, to the viscosity constant, and, for geometrically 
similar systems, to the linear dimensions of the body. To the domain of this "linear resistance" — which has aroused 
much interest, especially within recent years, on account of some important experimental applications — has to be 
opposed the limiting domain of comparatively large velocities, for which the so-called "velocity square law" holds 
with ver>' good approximation. In this latter domain, whicli embraces nearly all the important technical applications, 
the resistance is nearly independent of fluid viscosity, and is proportional to the fluid density, the square of the 
velocity, and — again for geometrically similar systems — to a surface dimension of the body. In this domain of the 
"square law" is included the important case of air resistance, because it is easy to verify, by the calculation of the 
largest density variations which can occur for the speeds we meet in aeronautics and airscrews, that the air compression 
can be neglected without any sensible error. The influence of the compression first becomes important for velocities 
of the order of the velocity of sound. In fact, experiments show that the air resistance, in a broad range from the small 
speeds at which the viscosity plays a role up to the large speeds comparable to the velocity of sound, is proportional 
to the square of the velocity with very good approximation.' In general, fluid resistance depends upon the form and 
the orientation of the body in such a complicated way that it is extraordinarily difficult to predetermine the flow to a 
degree sufficient for the evaluation of the resistance of a body of given form, by a process of pure calculation, as can be 
done by aid of the Stokes formula in the case of very slow motions. We also will not succeed in this paper in reaching 
such a solution, but wUl still make the attempt to give a general view of the mechanum of fluid Te>,isUince within the 
limiix of the square law. 

We can state the problem of fluid resistance in the following somewhat more exact way. 

Since the time of the fundamental considerations of Osborne Reynolds on the mechanical similitude of flow ' 
phenomena of incompressible viscous fluids of different density and viscosity and — under geometrical similitude — for 
different sizes of the system considered, it is known that the resistance phenomenon depends upon a single parameter 
which is a certain ratio of the above-mentioned quantities. Thus the fluid resistance of a body moving with the uniform ' 
velocity Uin an incompressible unlimited fluid may be expressed by a formula of the form * 

(I) w=nlUf(^^'J 
where 

It is the viscosity constant 
p the fluid density 

I a definite but arbitrarily chosen linear dimension of the body, and f^Eli^ a function of the single variable 
R—^^- We will call "Reynolds' parameter" the quantity R which has a zero dimension. 

Theorj' and experiment show that for very small values of R — that is, for low velocities, or small bodies, or great 
viscosity — the function / (7?) is very nearly constant; the resistance coefficient of the Stokes formula corresponds to the 
limiting case of f(R) for R=0. The square law corresponds to the limiting case of R=oa . We approach this latter case 
the more nearly the smaller the viscosity ix, so that in the limiting case of 2?=ec, the fluid can be considered as 
frictionless. And we can ask ourselves, to what limiting configuration does the flow of the viscotis fluid around a solid 
body tend when we pass to the limiting case of a perfect fluidf This is, according to our view, the fundamental point of 
the resistance problem. 

The fact that we obtain in this case a resistance nearly independent of the viscosity constant — since according 
to formula (/) this corresponds to the square law — allows us to conjecture that in this limiting case the resistance is 
determined by flow tj'pes such as can occur in a perfect fluid. 

' This noto is the translation of the paper of Th. v. K&rmin and II. Rubach published in " Physikalischn Zcitschrift," Jan. 15, 1912. The 
author has considered it necessary to add here this complete translation, on account of the importance of the new conceptions of Th. v. K&rm&n. 

> The actual noto constitutes a more complete ex-position of two notes of Th. v. KArm/in published in the "Nachrichten der Kgl. Gesellsch. 
der Wiss. tu Gottingen," contaming his hydrodynatnical researches upon the stability ol vorlex systems and the conclusions concerning fluid 
resistance obtained from the latter. The experiments here discussed and the measurements given have a provisional character; exact measurements 
are expected in connection with an intended dissertation of H. Rubach. 

> See, for example, the Interesting experiments of O. Fdppl on the validity of the square law for air resistance, "Dissertation," ;Aachen, p. 40 
(also " Jahrbuch der MotorluItschia"-Stadiengesellschatt, 1911). 

« Compare Lord Rayleigh, Phil. Mag., vol. 21, p. 70S, 1911. 
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It is now certain that neither the so-called "continuous" potential flow, nor the "discontinuous" potential flow 
discovered by Kirchhoff and v. Eelmholtz, can express properly this limiting case. Continuous potential flow does not 
cause any resistance in the case of uniform motion of a body, as may be shown directly by aid of the general momentum 
theorem; the theory of the discontinuous potential flow, which, in relation to the resistance problem has been dis- 
cussed principally by Lord Rayleigh,' leads to a resistance which is proportional to the square of the velocity; the cal- 
culated \ iilucB do not, however, agree with the observ ed ones. And, independent of the insufficient agreement between 
the numerical values, the hypothesis of the "dead water," which, according to this theory ought to mo^•e with the 
body, is in contradiction to nearly all observations. It ia easy to see by aid of the simj)lost experiments that the flow, 
when referred to a system of coordinates moving with the body, is not stationary, as assumed in this theorj'. Further- 
more, in the theory of discontinuous potential motion, the suction effect behind the body is totally missing, while in 
the dead water, which extends to infinity, we have everywhere the same pressure as in the undisturbed fluid at a great 
distance from the body. But according to recent measurements, in many cases the suction effect is of Urst importance 
for the resistance, and in any case contributes a sensible part of the last. 

The reason why in a perfect fluid the discontinuous potential flow, although hydrodynamically possible, is not realized 
is without any doubt the instability of the surfaces of discontinuity, as has already been recognized by v. Ilelmholtz 
and specially mentioned by Lord Kelvin.' A surface of discontinuity can l)e considered as a vortex sheet; and it can 
be shown in a quite general way that such a sheet is always unstable. This can also be obser\'ed directly; observation 
shows that vortex sheets have a tendency to roll themselves up; that is, we see the concentration around some points 
of the vortex intensity of the sheet originally between them. This observation leads to the question: Can there exist 
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stable arrangements of isolated vortex filaments, which can be considered as the final product of decomposed vortex 
sheets? This question forms the starting point of the following investigations; it will, in fact, appear that at least for 
the simplest case of uniplanar flow, to which we will limit ourselves, we will be led to a "flow picture" which in all 
respects corresponds quite well to reality. 

TSra INVESTIGATION 07 STABIUTV. 

We wiU investigate the question whether or not two parallel rows of rectilinear infinite vortices, of equal strength 
but of inverse senses, can be so arranged that the whole system, while maintaining an invariable configuration, will 
have a uniform translation and be stable at the same time. It is easy to see that there exist two kinds of arrangements 
for which two parallel vortex rows can move with a uniform and rectilinear velocity. The vortices may be placed 
one opposite the other (arrangement a, fig. I), or the vortices of one row may be placed oppo.site the middle points 
of the spacing of the vortices of the other row (arrangement b). In the case of equality of spacing of the vortices in 
both rows, as a consequence of symmetry for the two arrangements a and b, it appears that each vortex has the same 
velocity in the sense of the X axis, and that the velocity in the sense of the >"axis is equal to zero. We have to answer 
the question, which of these two arrangements is stable? 

To illustrate first by a simple example the method of the investigation of stability, we will start with the con- 
sideration of an infinite row of infinite vortices disposed at equal distances I and having the intensity f , and will study 

' On the resistance ol Bulds, Mathematical and Physical Papers, Vol. I, p. 287. 

I Mathematical and Physical Papers, Vol. rv, p. 21S. This paper contains a detailed critique of the theory of discontinuoas motion. 
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the stability of such a system. If we designate by ip, y^, the coordinates of the p—th vortex, and by x,, jy,, the 
coordinates of the q—th the velocity impressed on the latter vortex by the former is given by the formula 

u L Vp-yq 

"pq— ^ 



2t (Xp-x,)H(!?p-!/,)» 

These formula express tie fact that each vortex communicates to the other a velocity which is normal to the line 
joining them and is inversely proportional to their distance apart. Therefore the resultant velocity of the q—th vortex 
due to all the vortices is equal to 

00 

di 2;rZ^'(a-p-Z,)»+(Vp-J/,)' 

P— — QO 

00 

dt 2TZj(2p-i,)'+(!/p-y,)' 

p— — oo 

whore p=q is excluded from the summation. If now the vortices are disturbed from their equilibrium position, the 
small displacements being fp, ijp, the vortex velocities can be developed in terms of these quantities, and we will 
be brought to a system of differential equations for the disturbances {p, tj^,, i. e., for small oscillations of the system. 
Let us accordingly put 

!'p=''p 

and, neglecting the small quantities of higher orders, we will get 

00 

dt ~2tZj(P-9)'P 

B = — 00 

dt ~2x/ l^p-ol'P 

P— — 00 

The differential equations so obtained, which are infinite in number, are reduced to two equations by the sub. 
stitution 

These two equations are 

dt ~''°2ir/ I pH' 
dt ^°2ir^| p»P 

P— — OO 

with p^o 

The physical meaning of this substitution is easy to see: we consider a disturbance in which each vortex undergoes 
the same motion only with a different phase <p. Under such conditions we have to do with a wave disturbance and the 
system will be called stable, when for any value of <p, that is, for any phase difference between two consecutive vortices, 
the amplitude of the disturbance does not increa.«e \vith the time. 
Let us introduce the notation 

, s_iLV^«i?frJ- f Vl cos(py)-l 
*W-2^2_| P'P ~tP2j 

The foregoing equations then take the form 

Let us put Jo and ijo proportional to a*-' ; we will then find for each value of <p two values for X, that is 

X=±'c(*') 
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It follows that the vortex eystem considered is unstable for any periodic disturbance, because there is always 
present a positive real value of X, that is, the disturbance is of increasing amplitude. 

Appljdng this method in the case of two vortex rows we will find that the arrangement a, that is, the symmetrical 
arrangement, is likewise unstable, but that for the arrangement b there exists a value of the ratio hjl (h is the distance 
between the two rows, / is the distance between the vortices in tlie row) for wliich the system is stable. 

In both cases X can be brought to the form 

where A, B, Care functions of the phase difference v>. The system will be stable if (C^ — A') is positive for any val ue of 
For the symmetrical arrangement a, the functions A, B, Care expressed by the formiilx : 

" 00 

pH'—K' , \ ^ l—coep<p- 



1 pH^-K' , \^ /-coe 



V 

CO 



B(a,)= y^-P^—J^ sin (p<p) 



p-X 



But for (p=x we get 



C(.)=^[ctgh^ Cf )-t.h»(^)] 



so that this arrangement is unstable for any values of h and I. 
For the unsymmetrical arrangement 6 we find 



.1 . \^ (p+iyP-h' .\^ l-co8{p^) 

^(•P)- ^J[{j,+iyp+h^f ZJ p'P 



p-0 
p— 0 

We see now that C(ir)=o, so that in the place where (fi—ir, A must also be equal to zero, because, on account of tlie 
double sign, X takes a positive real value. This brings us to the condition 

[(p+iyP+Vf- (UpTlW 



But 



and 



p. 



p-0 



so that, as the necessary condition of stability we find the relation 
and for the ratio hjl we find the value 



coshy^=V2 



hll=0,283.... 



\ 
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For a certain value of the wave length of the disturbance, corresponding to »;=7r, we get X=o, that is, the system 
is in a neutral state. But it can be shown by calculation that our system is stable for all other disturbances. This 
unique disturbance has to be tested by further investigations. It can. however, be seen that a zero value for X must 
appear, because only one stable coufiguration exists. If this were not so, we would find for l/h a finite domain of sta- 
bility.' 

THE "flow picture. " 

The consideration of the question of stability has brought us to the result that there exists a particular config- 
uration of two vortex rows which is stable. The vortices of both rows have then such an arrangement that the 
vortices of one row are placed opposite the middle of the interval between the vortices of the other row, and the 
ratio of the distance h between the two rows to the distance I between the vortices of the same row has the value 



The whole system has the velocity 



which can also be written 



or, introducing th« value of hjl found by the stability investigation, we get 

The flow is given by the complex potential potential, flow function) 

if Bin(2„-2)j 

x=f+i>f'—2i Ig — — * 

sin («o+z)f 

where 

Zo-4+2 

By aid of this formula we have calculated the corresponding streamlines and have represented them in Fig. 2. 
We see that some of the streamlines are closed curN'ea around the vortices, while the others run between the vortices. 

On the other hand, we have tried to make %Tsible the flow picture behind 
a body, e. g., a flat plato or circular cylinder, moved through immobile 
water, by aid of lycopodium powder sifted on the surface of the water, and 
to fix these pictures photographically (exposure one-tenth of a second). 

The regularly alternated arrangement of the vortices can not be 
doubted. In most cases the vortex centers can also be well determined; 
sometimes the picture is disturbed by small "accidental vortices" pro- 
duced in all probability by small vibrations of the body, which in our pro- 
visional experiments could not be avoided. We had a narrow tank whose 
floor was formed by a band running on two rolls, and the bodies tested 
were simply put on the moving band and carried by it. It is to be expected 
that by aid of an arrangement especially made for the purpose much more r^:ular flow pictures could be obtained, 
while in the actual experiments the flow was disturbed on the one hand by the vibrations of the body and on the 
other by the water flow produced by the moving band itself. 

The alternated arrangement of the vortices rotating to the right and to the left can only be obtained when the 
vortices periodically run off first from one side of the body, then from the other, and so on, so that behind the body 
there appears a periodic motion, oscillating from one side to the other, but with such a regularity, however, that the 
frequency of this oscillation can be estimated ^^dth sufficient exactness. The periodic character of the motion in the 
so-called "vortex wake" has often been observed. Thus, Bernard ^ has remarked that the flow picture behind a 
narrow obstacle can be decomposed into vortex fields with alternated rotations. Also for the flow of water around 
balloon models the oscillation of the vortex field has been observed.' Finally, v. d. Borne * has observed and pho- 
tographed recently the alternated formation of vortices in the case of air flowing around different obstacles. The 

I From a mathematical standpoint our stability investigation may be considered as a direct application of tlie theorems of Mr. O. Topliti on 
Cyclanten with an infinite number of elements, wliich be has in part published in two papers (Gottingeo Naehrichten, 1907, p. 110; Math. Auualen 
1911. p. 361), and in part been so kind as to conununicate personally to us. 

> Comptes Rcndus, Paris. 14S, S.W, 190S. 

» Technical report of the Advisory Commlltoo for Aeronautics (British), WlO-il. 

' Undertaken on the initiative of the representatives of aeronautical science in Gottingon, November, 1911. 
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phenomenon could not be explained until now; according to our stability investigation the periodic variations appear 
as a natural ronserjuence of the insfcibility of tho symmetrical flow.' 

It is also very interesting to observe how the stable configuration is established. When, for example, a body is 
set in motion from rest (or conversely, the stream is directed onto the body) some kind of "separation layer" is first 
formed, wliich gradually rolls itself up, at first symmetrically on both sides of the body, till some small disturbance 
destroys the symmetry, after which the periodic motion starts. The oscillatory motion is then maintained corresponding 
to the regular formation of left hand and right hand vortices. 

We have also made a second series of photographs for the case of a body placed at rest in a imiform stream of 
water. For this case the flow picture can be obtained from Fig. 2 by the superposition of a uniform horizontal velocity. 
We will then see on the lines drawn through the vortex centers perpendicular to the stream direction, some ebbing 
point where the stream lines intersect and the velocity is equal to zero. However, in the same way as the motion is 
affected by the \-ibrations of the experimental body in the case of the motion of a body in the fluid, so in this case the 
turbulence of the water stream gives rise to disturbances. 

As to the quantitative agreement attained by the theory, it must be noted that our stability conditions refer to 
infinite vortex rows, so that an agreement of the ratio /(// with the measured values is to be expected only at a certain 
distance from the body. The measurements on the photographs show that the distance I between vortices in a row is 
very regular, so that I may be measured satisfactorily, but per contra the distance h is much more variable, because 
the disturbance of the vortices takes place principally in the direction normal to the rows, that is, the latter undergo 
in the main .transverse oscillations. The best way to determine the mean positions of the centers of the vortices would 
be by aid of cinematography, but we can also, without any special difficulty, find by comparison the mean direction 
of each vortex row directly from photographs. So in the case of the photograph of a circular cylinder 1.5 cm. in diameter, 
when making measurements beyond the first two or three vortex pairs we have found the following mean values for 
h and I 

A=1.8 cm.; 1=6.4 cm. 

So that for the ratio h/l we obtain the value 

hll=0.28. 

For the flow around a plate of 1.75 cm. breadth we found 

/(=3 cm.: /=9.8 cm. 

.Accordingly 

A/i=0.306. 

The agreement with the theoretical value 0.283 is entirely satisfactory. 

For the first vortex pair behind the body, h/l comes out sensibly larger, somewhere near A/i=0.35. But in the first 
investigation of KArmdn, mentioned at the beginning of this paper, the stability of the vortex system was investigated 
in such a way that all the vortices with the exception of one pair were maintained at rest and the free vortex pair con- 
sidered oscillating in the velocity field of tlie others. Under such assumptions it was found that A/'=1/t arc cosh V3^= 
0.36. We therefore think that the conclusion can be drawn, that in the neighborhood of the body, where the vortices 
are even more limited in their displacements, the ratio h/l is greater than 0.283 and approaches rather the value of 0.36. 

APPLICATION OP THE MOME.VTUM THEOREM TO THE CALCULATION OP FLUID RESISTANCE. 

Let us assume that at a certain distance behind the body there exists a flow differing but slightly from the one 
of stable configuration which we have established theoretically in the foregoing, but that at a distance in front of the 
body, which is great in comparison with the size of the body, the fluid is at rest — as it is quite natural to assume. 
We will then be brought by the application of the momentum theorem to a quite definite expression for the resistance 
which a body moNing with a uniform velocity in a fluid must experience. Practically, by such a calculation for the 
uniplanar problem, we will obtain the resistance of a unit of length of an infinite body placed normally to the plane 
of the flow. 

We will use a system of coordinates moving with the same speed u as the vortex system behind the body. In 
this coordinate system, according to our assumptions, at a sufficient distance from the body the vortex motion behind 
the body as well as the fluid state in front of the body will be steady, and we will have, when referred to this system 
of Coordinates, a uniform flow of speed -u in front of the body, but behind the body the velocity components will 
be expressed by 

-w+5^and-5^ 

where ^ is the real part of the complex potential 



if sin (zo+z) j 
Bin (zo-z) J 



> The tone that Is emitted by a stick rapidly displaced in air is Hxcd by this periodicity, to which Prof. C. ICunge has already drawn our 
attention. 
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The body itself has, relative to this system of coordinates, the velocity U—u. where U is tlie absolute vcU)city of the 
body. If we designate by I the distance between the vortices of one row, there must take place, as a consequence 
of the displacement of the body, in the time T=ll{ U—u), the formation of a vortex on each side of the body. We 
will calculate the increment of the momentum, along the X axis, in this time interval 7' (that is, between two instants 
of lime of identical flow state) atid for a part of the flow plane, which we define in the following way (see fig. 3). On 
the sides the i>Iane [xirtion considered is limited by the two parallel straight lines j/=±i).' in front and behind, by 
tw<i straight lines 2=f'on8t. disposed at distances from the body which are great in comparwon with the size of the 
body, the line behind the body being drawn so as to pass through the |X)int half way between two vortices having 
inverse rotation. When the boundary lines are sufliciently far from the body we can consider the fluid velocitieB at 
those lines as having the values indicated in the foregoing. 

For a space with the boundaries indicated above the relation must exist that the momentum imparted to the 

body ^ gWdl (where Wis the resultant fluid resistance) is equal to the difference between the momentum contained 

in the space considered at the times t=r and t=T4-7'and the sum of the inflow momentum and the time integral 
of the pressure along the boundary lines. If we thus consider as exterior forces the force — W and the pressure, which 
act on the whole system of (luid and solid, they must then correspond to the increment of the momentum — that is, 
to the excess of momentum after the time 7' less the inflow momentum. 
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We will calculate these momentum parts separately. The excess of momentum after the time T is equal to the 
difference of the values that the double integral pjju (x, y) dx dy takes at the times t=T and «=r4- T. But the 
time interval has been chosen in such a way that the state of flow is identical, Mrith the difference that the body has 
been displaced through the distance /=( U—u) T. The double integral reduces thus to the difference of the integrals 
taken over the strips ABCD and A'B'C'D' both of breadth I. For the strip A'B'C'iy the fluid speed can be taken 

equal to — u for the strip ABCD equal to— go that we get 



■■■■//: 



^dxdy 



If we pass to side boundaries having if=aB , we obtain for 7, the very simple expression 

which can also be obtained directly by the application of the general momentum theorem to vortex systems. 

We will unite in one single term the inflow momentum and the time integral of the pressure, because in such a 
way we will be led to more simple results. If we consider a uniplanar steady fluid motion with the velocity com- 
ponents w {x, y) and v (i, y) and consider a fixed contour in the plane, the inflow momentum in a unit of time in the 
direction of is expressed by the closed integral pj{t?dy—uvdx) where u, v are the velocities on the contour. The 
pressure gives the resultant jpdy along the X axis, but since for a steady flow the relation 

u'+v' 

p=ConBt— p — 2 — 

must hold, we thus obtain for the sum of both integrals, multiplied by T 



I.,= Tjf,(u^dy-uvdx)+Tfp dy 
= 7'p J ( ^2^ ~ uvdx^ 
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Or, introducing the complex quantity, 
we get 



w=u—iv= ST — r"=*\=5C~ 



where Im is to be understood as the complex part of the integral. 
If we put for the contour 

u=—u-\-u' 

then the terms in will at once bo eliminated, and also the terms in u on account of the e((uality of the inflow and 
outflow; and there will remain only the terms in u'^ and u'v'. The latter will give a finite value only for the boundary 
line passing through the vortex system (AD in fig. 3). Passing to 17=00 , we get 



l2=TpIm 



and integrating along AD we get 



But 



^Tplm 



r rm 



if COB 



2irx 



I cosh 



80 that, integrating and introducing the values 



x(-iw)=- 



l.-Tp^-l — 



where u again has been written for^ ^S^^y ' 

Thus the total momentum imparted to the body is 



If for the mean value of '/t^Wdt we write W (as the time mean value of the resistance) we will obtain with 
T=l/( U-u) the final formula 



(ID 



W=pt-j(U-2u)+p-l^^ 



The fluid resistance appears here expressed by the three characteristic constants f, h, I of the vortex configuration 
(aa u is expressed by the last). In the deduction of this last formula we did not take account of the stability condi- 
tions, BO that this formula applies to any value of the ratio h/l. If we assume the vortices in the row to be brought 
all close together so that they are uniformly distributed along the row, but in such a way that the vortex intensity per 
unit of length remains finite, we thus pass to the case of continuous vortex sheets. In this case f/i= U, but iyi=0 



and " =-^' so that the fluid resistance disappears. 



The discontinuous potential flow of v. Helmholtz thus does not 
give any resistance when the depth of the dead water remains tinite, as can also be shown from general theorems. 

THE FORMULAE FOR nVTD RESISTANCE. 



Let us now apply to our special case the general formula we have just found, introducing the rdilations between 
f and u, and h and I according to the stability conditions. For the speed u we have ^ 
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further. 

A//=0,283 

80 that we get 



If we introduce, as is ordinarily done, the resistance coefficient according to the formula 

IF=^„ pdW 

where rf is a chosen characteristic dimension of the body, to which we refer the resistance, we will obtain ^„ expressed 
by the two ratios u/ U and Ijd in the following way 



(III) *„=[o,799-^-0,323(-^)'] 



We have thus obtained the resistance coefficient — which before could be observed only by resistance measurements — 
expressed by two quantities which can be taken directly from the flow phenomenon, viz, the ratio 

u _Velocity of the vortex s ystem 
'D~ Velocity of the body 

and 

J _ Distances apart of the vortices in one row 
d~ Reference dimension of the body 

Both quantities, corresponding to the similitude of the phenomenon, within the limits of validity of the square law 
can depend only upon the dimension of the body. 

The.'^e two quantities can be observed very easily experimentally. The ratio Ijd can be taken directly from photo- 
graphs, while the ratio uj U can be found easily by countinR the number of vortices formed. If we designate by T the 
time between two identical flow sUites we can then introduce the quantity lo= OT, which is the distance the body 
moves in the period T. This quantity must be independent of velocity for the same body, and the ratio Ijlg for similar 
bodies must also be independent of the dimensions of the body but determined by the shape of the body. Remember- 
ing that T=ll( U—u), we then find between «./ f/and l/lo the simple relation 

lo 

By some provisional measurements we have proved the similitude rule and afterwards calculated the resistance 
coeflScient for a flat plate and a cylinder disposed normal to the stream, for the purpose of seeing if the calculated values 
agreed with the air resistance measurements, at least in order of magnitude. 

Our measurements were made first on two plates of width 1.75 and 2.70 cm. and 25 cm. length, and we have meas- 
ured the period T and calculated the quantity lo=UT for two different velocities. We have used a chronograph 
for time measurements and the period was observed for each vortex row independently. Thus was found for the 
narrower plate 

U=W.O cmfsK 15.1 m/««c 
T=1.26 sec 0.805 sec 



iot the wider plate 



f/r=12.6 cm 12.1 cm 

U=9.G cmjsec 15.5 cmjsec 

r=1.99»ec. 1.20 sec. 

UT=l9.l sec. 18.6 sec. 

Mean value i7r=18.8 cm 



The ratio of the plate width is equal to 

and the ratio of the quantities lo= UT \s equal to 



j-y3=1.54 



^-5J-l 52 
12. 3-^-^^ 



So that the similitude rule is in any case confirmed. 

A circular cyJiadeLoLl.5 em. diameter was also teaated at two speeds. We found the values 

{7=11.0 cm/8«c Xb.Z cmjsec 
7'=0.6e sec. 0.48 sec. 

UT=T.Z cm 7.5 cm 



Mean value UT=7.4 cm 
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Knowing the values of lo= UT we can calculate for the plate and the cylinder the speed ratio u/Z7. Thus, 

for the plate u/ f7=0.20. 
for the cylinder w/ (7=0.14 

and with the values of { indicated before wo have 

for the plate J/rf=5.5 
for the cylinder Z/rf=4.3 

where d is the plate width or cylinder diameter. We thus find the resistance coefficients 

for the plate <p^.=0.80 
for the cylinder ^,=0.4G 

The resistance measurements of Foppl' have given for a plate with an aspect ratio of 10:1 the resistance coefficient 
^^=0.72 and the Eiffel- measurements, for an aspect ratio of .50:1 ; that is, for a nearly plane flow, the value ^„=0.78. 
Further, Foppl has found fora long- circular cylinder if/ „=0.4o, so that the agreement between the calculated and meae- 
lured re.sistance coefficients mu.st be considered as fully satisfactorj-. 

The theoretical investigations here developed ought to be extended and completed in two directions. First, 
we have limited ourselves to the uniplanar problem; that is, to the limiting case of a body of great length in the direc- 
tion normal to the flow. It is to be expected that by the investigation of stable vortex configumtions in space we will 
al« be brought to a better understanding of the mechanism of fluid resistance. However, the problem is rendered 
difficult by the fact that the translation velocity of cur\-ed vortex filaments is not any longer independent of the size 
of the vortex section, because to an infinitely thin filament would correspond an infinitely great velocity. Never- 
thele.S3, it must not be considered that the extension of the theory to the case of space would bring unsurmountable 
difficulties. 

Much more difficult appears the extenaidn of the theory in another direction, which really would first lead to a 
complete understanding of the theory of fluid resistance, namely, the e\aluation by pure calculation of the ratios 
•'/rfaiuj u/U, which we have found from flow observations, and which determine the fluid resistance. This problem 
can not be solved without investigation of the process of vortex formation. An apparent contradiction is brought out 
ky th ! fact that we have .-sod only the theorems established for perfect fluids, which in such a fluid (frictionless fluid) 
iJrj vortices can be formeJ. This contradiction is explained by the fact that we can everywhere neglect friction except 
at-, \e •;uriace of the lioiiy. It con be shown that the friction forces tend to zero when the friction coefficient decreases, 
but the vortex intensity remains flnite. If we thus consider the perfect fluid as the limiting case of a viscous fluid, 
Jiui the law of vortex formation mut' be limited by the condition that only those fluid particles can receive rotation 
which have been in contact with the surlace of the body. 

This idea appear first, in a perfectly clear way, in the Prandtl theory of fluids having small friction. The Prandtl 
theory investigates those phenomena which take place in a layer at the surface of the body, and the way in which the 
separation of the flow from the surface of the body occurs. It we could succeed in bringing into relation these inves- 
tigations on the method of separation of the stream from the wall with the calculation of stable configuration of vortex 
films formed in any way whatever, as has been explained in the foregoing pages, then this would evidently mean 
great progress. Whether or not this would meet with great difficulties can not at the present time be stated. 

' S«e the work ot O. Foppl already mentioned. 

> G. Eiflel, "La Resistance de l'.\ir et l'.\viation," p. 47, Paris, IVIO. 
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